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PREFACE 

An  important  consideration  for  the  chemical  engineer  is  his 
ability  to  predict  reliable  values  for  various  properties  of  chemical 
substances  that  are  needed  for  design  equations.   Because  of  the 
large  number  of  combinations  of  substances,  compositions,  pressures, 
temperature,  etc., that  may  be  encountered,  it  is  often  neither 
feasible  nor  desirable  to  make  experimental  measurements  of  such 
properties.   The  need  for  reliable  correlations  has  become  more 
acute  with  the  increasing  use  of  high-speed  electronic  computers  in 
designing  chemical  plants.   Such  design  requires  suitable  mathe- 
matical expressions  for  properties  as  a  function  of  operating 
conditions.   Even  when  experimental  data  in  tabulated  form  are 
available  it  is  most  conveniently  introduced  into  the  computer  in 
the  form  of  soundly  based  theoretical  equations. 

In  the  long  term  the  only  satisfactory  approach  is  one  that 
is  firmly  based  on  molecular  considerations.   Theories  for  dilute 
fluids  are  well  established.   The  kinetic  theory  of  gases  is  capable 
of  describing  the  thermodynamic  and  transport  properties  accurately. 
However,  for  dense  fluids  there  is  still  no  theory  which  can  parallel 
the  success  of  the  dilute  gas  theory.  Although  there  are  several 
rigorous  molecular  theories  for  liquids,  none  have  yet  been  developed 
to  a  stage  that  would  yield  numerical  results. 

Theoretical  study  of  transport  properties  is  more  difficult 
than  that  of  thermodynamic  properties.   In  thermodynamics,  all 
equilibrium  properties  can  be  obtained  once  the  partition  function 
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of  the  system  is  kno^^m.   However,  for  transport  processes  each  of 
the  transport  properties  has  to  be  formulated  separately.   This 
situation  occurs  because  the  thermodynamic  equilibrium  state  is 
unique,  whereas  there  are  many  types  of  non-equilibriian  state. 

This  dissertation  considers  corresponding  states  correla- 
tions of  transport  properties  of  fluids,  particularly  liquids.   This 
principle  has  previously  proved  of  great  value  to  engineers  in 
predicting  thermodynamic  properties.   Chapter  1  provides  a  brief 
survey  of  the  present  status  of  theories  of  transport  properties  of 
dense  fluids.   Chapters  2-4  contain  a  detailed  study  of  the  free 
volxome  theory  of  viscosity  and  self-diffusion  coefficient,  with 
particular  emphasis  on  developing  corresponding  states  relations 
for  the  parameters  involved.   The  next  two  chapters  discuss  a  more 
direct  corresponding  states  treatment,  and  a  new  theory  is  proposed 
for  polyatomic  molecules. 

The  author  is  greatly  indebted  to  Dr.  K.  E.  Gubbins,  chairman 
of  his  supervisory  committee,  for  his  interest,  stimulation,  advice 
and  encouragement  during  the  course  of  this  research.   He  is  also 
grateful  to  Dr.  T.  M.  Reed  for  helpful  discussions.  Finally,  he 
expresses  his  sincere  appreciation  to  Dr.  T.  M.  Reed,  Professor 
R.  D.  Walker  and  Dr.  R.  G.  Blake  for  serving  on  the  committee. 
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CORRESPONDING  STATES  RELATIONSHIPS  FOR  TRANSPORT 
PROPERTIES  OF  PURE  DENSE  FLUIDS 

By 

Min  Jack  Tham 

June,  1968 

Chairman:   K.  E.  Gubbins 

Major  Department:   Chemical  Engineering 

Two  approaches  have  been  considered  for  developing  corres- 
ponding states  correlations  for  dense  fluid  transport  properties. 
The  first  is  based  on  the  free  volume  theory  for  viscosity  and  self- 
diffusion  coefficient.   This  theory  has  been  improved  to  account 
quantitatively  for  the  nonlinearity  of  constant  volume  plots  of 
ln(T]/T   )  vs.  T~   and  InCD/T''"'  )  vs.  t"''"  for  simple  liquids,  by  the 
introduction  of  the  temperature  dependence  of  V  and  the  volume 
dependence  of  E  .   Corresponding  states  relations  are  presented  for 
these  two  quantities.   The  improved  free  volume  theory  accurately 
describes  the  viscosity  and  self-diffusion  coefficient  behavior 
of  liquids  composed  of  simple,  nonpolar  molecules  that  may  be  expected 
to  approximately  obey  a  Lennard- Jones  [6,12]  potential  law,  for  the 
density  range  of  p  >  2p  .    It  is  in  this  range  that  most  other 
theories  fail.   The  only  parameters  required  for  the  viscosity  and 
self-diffusion  coefficient  prediction  are  melting  and  critical  para- 
meters, together  with  one  experimental  value  of  viscosity  (or  self- 
dif fusivity) . 

XX 


Direct  corresponding  states  correlations  of  transport  pro- 
perties have  also  been  studied,  and  provide  a  more  general  if  slightly 
less  accurate  method  of  prediction.   The  simple  corresponding  states 
principle  derived  by  Helfand  and  R.ice  from  the  time  correlation  func- 
tion expressions  of  statistical  mechanics  has  been  shown  to  apply 
closely  for  the  transport  properties  of  monatomic  molecules  over  the 
entire  range  of  temperature  and  density  conditions  for  which  measure- 
ments are  available.   In  this  simple  form,  however,  it  was  found  not 
to  apply  to  polyatomic  fluids  in  the  dense  phase.   After  careful  re- 
examination of  the  four  assumptions  on  which  the  derivation  of  the 
simple  correspondence  principle  is  based,  it  is  shown  that  hindered 
rotation  at  high  density  is  the  most  probable  cause  for  the  deviation 
of  polyatomic  fluids  from  the  correspondence  principle. 

A  correspondence  principle  for  polyatomic  molecules  has  been 
derived,  taking  into  account  the  effects  of  hindered  rotation  and 
making  use  of  a  simple  harmonic  oscillator  model.   The  proposed 
equations  have  been  tested  exhaustively  for  the  viscosity,  thermal 
conductivity  and  self-diffusion  coefficient  of  a  large  number  of 
polyatomic  fluids  over  a  wide  range  of  temperature  and  density  condi- 
tions.  The  molecules  tested  include  saturated  hydrocarbons  from 
methane  to  nonane,  spherical  molecules  such  as  neopentane  and  carbon 
tetrachloride,  flat  molecules  such  as  benzene,  and  simple  diatomic 
molecules  such  as  nitrogen.   The  improved  correspondence  principle 
predicts  the  transport  properties  of  saturated  liquids  within  experi- 
mental error.  For  high  pressures  the  average  percentage  errors  of 
the  predicted  viscosity  and  thermal  conductivity  are  slightly  higher. 

XX  i 


It  is  also  shown  that  the  proposed  model  can  account  in  a  reasonable 
way  for  thermodynamic  properties  of  polyatomic  fluids. 
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CHAPTER  1 
INTRODUCTION 

In  this  chapter  a  brief  survey  and  evaluation  is  presented  of 
the  principal  theories  and  empirical  correlations  that  have  been  pro- 
posed to  explain  dense  fluid  transport  properties.   For  pure  fluids 
all  transport  properties  are  functions  of  temperature  and  density, 
namely,  r]  =  rj(T,p),  /c  =  /c(T,p),  D  =  D(T,p)  and  >,  =  \(T,p).      None  of 
the  present  theories  provide  a  satisfactory  means  of  predicting 
transport  properties  over  the  entire  range  of  temperature  and  density 
conditions  for  even  quite  simple  fluids.   Also  assumptions  of  dubious 
validity  must  be  made  in  most  cases  in  order  to  be  able  to  make  any 
comparison  with  experiment  at  all.   In  the  following  survey  the  more 
rigorous  approaches  are  first  described,  followed  by  simple  model 
theories  and  ending  with  the  empirical  correlations. 
1 . 1 The  Time  Correlation  Function  Theory 

In  the  last  decade  a  new  theory,  called  the  time  correlation 
function  theory  (also  known  as  the  fluctuation-dissipation  theory) 
has  been  developed  to  treat  transport  processes.   This  theory  treats 
transport  processes  in  terms  of  equilibrium  ensembles  whose  properties 
are  known.   In  a  sense  the  time  correlation  function  plays  a  similar 
role  to  the  partition  function  in  statistical  thermodynamics.   In 
statistical  thermodynamics,  all  thermodynamic  properties  of  any  parti- 
cular system  can  be  evaluated  if  the  partition  function  is  known.   In 
the  same  manner  the  transport  coefficients  can  be  obtained  if  the 
appropriate  time  correlation  functions  can  be  calculated.   However,  in 


one  respect  the  analogy  breaks  down.   In  statistical  thermodynamics  a 
single  partition  function  determines  all  the  thermodynamic  properties, 
whereas  in  transport  processes  different  time  correlation  functions 
are  needed  for  different  transport  processes.   A  good  review  of  the 
time  correlation  function  approach  to  transport  processes  is  given  by 
Zwanzig  ( 1) . 

The  time  correlation  function  expressions  for  the  transport 
coefficients  may  be  derived  by  several  methods.   A  simple  derivation 
involves  starting  with  expressions  analogous  to  the  Einstein  equation 
for  the  self-diffusion  of  Brownian  particle  (2).   The  following 
expressions  for  transport  coefficients  are  obtained. 
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For  dilute  gases  the  time  correlation  functions  may  be  evaluated  and 
the  expressioiB  for  the  transport  coefficients  can  be  reduced  to  the 
forms  obtained  by  solving  the  Boltzmann  equation.   However,  for  dense 
fluids  no  satisfactory  expressions  have  as  yet  been  obtained,  although 
several  approximate  equations  have  been  suggested  (3) .   One  approach 
has  been  to  assume  that  the  time  correlation  function  decays  exponen- 
tially.  Alternatively  one  may  perform  molecular  dynamics  calculations 
(4-6) .   In  this  method  a  very  large  amount  of  computation  is  involved; 
the  principal  interest  of  such  an  approach  is  to  provide  data  with 
which  to  compare  various  theoretical  expressions  for  the  time  correla- 
tion functions. 

The  time  correlation  theory  gives  a  formal  description  of 
transport  processes;  however  it  is  still  not  developed  to  a  sufficiently 
advanced  stage  to  be  of  practical  value  in  the  calculation  of  trans- 
port coefficients.   To  the  chemical  engineer,  the  most  interesting 
feature  of  the  theory  in  its  present  state  of  development  is  that  it 
provides  the  basis  for  a  rigorous  development  of  the  corresponding 
states  principle  for  dense  fluid  transport  properties.   This  is  dis- 
cussed in  a  later  section  of  this  chapter. 
1 .  2 The  Kinetic  Theory  of  Liquids 

The  kinetic  theory  of  dilute  gases  is  already  very  well 
developed  (7),  and  provides  a  method  of  accurately  estimating  their 
transport  properties.   It  is  therefore  reasonable  to  attempt  a  similar 
development  for  dense  gases  and  liquids.   Such  an  attempt  was  first 
made  by  Kirkwood  (8,9),  whose  approach  was  based  on  casting  the 


Liouville  equation  for  the  distribution  function  into  an  equation  of 
the  Fokker-Planck  type.   The  main  disadvantage  in  this  theory  is 
that  a  large  number  of  approximations  have  to  be  made  in  obtaining  the 
Fokker-Planck  equations  for  f    and  f    ,  the  time  smoothed  singlet 
and  pair  distribution  functions  respectively.   In  addition  to  these 
deficiencies,  Kirkwood's  method  has  drawn  much  criticism  for  assuming 
the  mean  momentum  change  for  a  collision  to  be  small.   The  equations 
derived  by  the  method  of  Irving  and  Kirkwood  (9)  for  transport  coeffi- 
cients of  liquids  are  as  follows  (10) 
Shear  viscosity       ,  „  oo 

T]  =^-^^    r'  ^  Z^(r)^,(r)ar  (1.8) 

^    15kT   ^     ar    o    z 

Bulk  viscosity      r   ^    °° 

Thermal  conductivity  „       . 

5k  niT        nVTTkT  °° 
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O      00 
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kT 
Self-diffusion  coefficient    D  =  — r  (1.11) 

where  °° 

C   =^p   rr^g(r)v^u(r)dr  (1.12) 

^o 
and     ■^o'^2  ~  ^^^  radial  functions  in  the  spherical  harmonic 
expansion  of  the  nonequilibrium  distribution 
function  g(r) 


A  modification  of  the  Kirkwood  theory  due  to  Rice  and  Allnatt 
(11,12)  attempts  to  avoid  the  assumption  that  the  mean  collisional 


momentum  change  is  small.   In  the  Rice-Allnatt  theory,  a  potential 
energy  function  is  assumed  in  which  the  molecule  has  a  hard  core 
together  with  an  outer  shell  which  interacts  with  a  soft  attractive 
potential.   A  hard  core  collision  leads  to  a  large  momentum  transfer 
and  tends  to  vitiate  the  Brownian  motion  approximation  which  is  inherited 
in  the  Fokker-Planck  equation.   These  authors  therefore  treated  the 
rate  of  change  of  the  distribution  function  f  due  to  hard  core  colli- 
sions by  means  of  an  Enskog-type  collision  term,  and  treated  the  rate 
of  change  of  the  distribution  function  due  to  motion  in  the  attractive 
field  of  surrounding  molecules  betv^7een  collision  by  Kirkwood's  approach. 

In  order  to  calculate  the  transport  coefficients  from  the 
Kirkx^70od  or  Rice-Allnatt  theory  one  must  have  accurate  values  for  the 
intermolecular  potential,  and  also  the  equilibrium  radial  distribution 
function  has  to  be  known  quantitatively.   The  latter  requirement  is 
especially  difficult  to  meet  since  the  experimental  radial  distribution 
function  has  not  been  determined  accurately. 

Making  use  of  the  Kirkwood  theory,  Zwanzig,  Kirkwood,  Stripp 
and  Oppenheim  (13)  were  able  to  calculate  the  shear  viscosity  and  bulk 
viscosity  for  liquid  argon  near  the  boiling  point.   Their  calculated 
value  for  shear  viscosity  was  r]  =  0.73x10   poise,  while  the  experimen- 
tal value  of  shear  viscosity  for  argon  at  the  same  temperature  is 

_3 
T)  =  2.39x10   poise.   Thus  the  calculated  value  is  in  error  by  roughly 

a  factor  of  three.   For  bulk  viscosity,  their  computation  yielded 

k;  =  0.36x10   poise.   However,  according  to  Naugle's  (14)  measurement 

of  bulk  viscosity  for  liquid  argon  at  approximately  the  same  temperature, 


-3 
K   =    1.7x10    poise.   The  discrepancy  between  calculated  and  experimen- 
tal values  is  more  than  a  factor  of  four.   Because  of  the  non- trivial 
mathematics  involved,  the  Kirkwood  theory  has  not  been  tested  for  the 
temperature  and  density  dependence  even  for  very  simple  molecules  such 
as  argon. 

The  Rice-Allnatt  theory  has  been  tested  by  Lowry,  Rice  and 
Gray  (15)  for  viscosity.   They  found  quite  good  agreement  between 
calculated  and  experimental  values  of  argon  at  densities  slightly  less 
than  the  normal  liquid  density  (Table  1.1).   To  compare  theory  and 
experiment  they  assiomed  u(r)  to  be  given  by  the  Lennard- Jones  [6,12] 
potential.   The  radial  distribution  functions,  g(r)  for  different 
temperatures  were  those  obtained  theoretically  by  Kirkwood,  Lewinson 
and  Alder  (16) .   The  radial  distribution  functions  g(r)  and  the 
derivative  of  u(r)  were  modified  by  introducing  an  empirical  parameter 
c^ ,  which  was  obtained  from  equation  of  state  data. 

TABLE  1.1 

TEST  OF   RICE-ALLNATT   THEORY  FOR 
SHEAR  VISCOSITY  OF   ARGON 


-3 

Density,    gem 

Temperature,   K 


Pressure,  atm. 


^ 


rj(calc) 
T](expt) 
7o  error 


r)xlO- 

,  poise 

1.12 

1.375 

128 

133.5 

185.5 

90 

50 

100 

500 

1.3 

0.9819 

0.9827 

0.9887 

0.9705 

0.727 

0.730 

0.771 

1.74 

0.835 

0.843 

0.869 

2.39 

13.0 

13.4 

11.3 

27.2 

Ikenberry  and  Rice  (17)  tested  the  Rice-Allnatt  theory  on 
the  thermal  conductivity  of  argon  for  the  same  temperature  and  density 
conditions  as  the  viscosity  values  shown  in  Table  1.1.   The  calculated 
values  were  not  entirely  theoretical  since  they  made  use  of  experi- 
mental friction  constants  calculated  from  the  self-diffusion  coeffi- 
cient data  of  Naghizadeh  and  Rice  (18) .   Close  agreement  was  found 
bet\-;een  theory  and  experiment  for  the  entire  temperature  and  density 
range  considered  as  shown  in  Table  1.2. 

TABLE  1.2 

TEST  OF  RICE-ALLNATT  THEORY  FOR 
THERMAL  CONDUCTIVITY  OF  ARGON 

XxlO^,  cal.cmT-'-secT-'-^c" 


_3 

Density,  gem. 

1.12 

1.375 

Temperature,  °K 

128 

133.5 

185.5 

90 

Pressure,  atm. 

50 

100 

500 

1.3 

X(calc) 

2.83 

1.85 

1.77 

1.84 

X(expt) 

2.96 

1.89 

1.86 

1.87 

%   error 

4.4 

2.1 

4.8 

1.6 

Davis,  Rice  and  Sengers  (19)  further  modified  the  Rice-Allnatt 
theory  by  using  the  square-well  potential  energy  function.   Because  of 
the  simple  form  of  this  function  the  Rice-Allnatt  theory  is  greatly 
simplified  and  the  evaluation  of  the  pair  correlation  function  reduces 
to  the  determination  of  g(o^)  and  zi^ '^  •      ^^i  order  to  calculate  the 
transport  coefficients  from  this  theory,  one  must  know  the  values  of 
the  parameters  O,,  a„,  £,  g(ci),  g(cr  )  and  have  equation  of  state  data. 


In  general  a  ,    a„  and  e   for  the  square-well  potential  are  available 
(20).   In  principle  g(a  )  and  g(a  )  can  be  calculated  from  equilibrivim 
statistical  mechanics.   However  there  is  still  no  satisfactory  method 
of  making  such  theoretical  calculations  accurately. 

Davis  and  Luks  (21)  estimated  the  values  of  g(0'^)  and  g(0'  ) 
from  equation  of  state  and  thermal  conductivity  data  of  argon  by 
solving  simultaneously  the  modified  Rice-Allnatt  expressions.   They 
then  calculated  the  viscosity  coefficient  and  self-diffusion  coeffi- 
cient of  argon  at  various  temperatures  and  pressures.   The  results  are 
shown  in  Figures  1.1  and  1.2.   The  discrepancy  between  theory  and 
experiment  for  viscosity  is  as  much  as  19%  and  for  the  self-diffusion 
coefficient  is  up  to  34%.   The  experimental  self-dif fusivity  data 
used  by  Davis  et  al.  as  shown  in  Figure  1.2  are  apparently  extrapolated 
values  obtained  from  the  self-diffusion  coefficient  data  of  Naghizadeh 
and  Rice  (18).   Luks,  Miller  and  Davis  (22)  have  also  made  use  of  the 
modified  Rice-Allnatt  theory  to  calculate  transport  coefficients  for 
argon,  krypton  and  xenon.   Discrepancies  between  theory  and  experiment 
become  larger  as  the  density  increases. 
1.3     The  Enskos  Theory  (7) 

Enskog  extended  the  dilute  gas  kinetic  theory  to  dense  fluids 
composed  of  hard  spheres.  As  a  gas  is  compressed  the  mean  free  path, 
which  for  dilute  gases  is  infinitely  large  with  respect  to  the  molecular 
diameter,  gradually  decreases  until  at  high  density  it  is  of  the  same 
order  of  magnitude  as  the  diameter  of  a  molecule.  During  this  process 
the  intermolecular  collision  frequency  tends  to  get  larger  due  to  the 
decrease  in  the  mean  free  path;  at  the  same  time  it  decreases  because 
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at  close  separations  the  molecules  tend  to  shield  one  another  from 
collisions  with  more  distant  neighbors.   The  net  change  in  collision 

frequency  was  found  to  differ  from  that  of  a  dilute  gas  by  a  factor 
^,  which  is  given  by 

X  =  1.0  +  0.625^-1)+  0.2869  (^  )  +0.115i;:^  )      (1.13) 
+  .  .  . 
where  b  =  —  ttNCt  =  the  rigid  sphere  second  virial  coefficient. 
Assuming  the  collisional  transfer  of '  momentimi  and  energy  between 
colliding  molecules  to  be  instantaneous,  Enskog  modified  the  Boltzmann 
equation  for  higher  densities  and  solved  it  to  obtain  the  following 
equations  for  transport  coefficients. 

b 


Self-diffusion  coefficient 


oD 


(pD)^   Vy 


(1.14) 


Shear  viscosity 


_o 
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+  0.8   +  0.761y 


(1.15) 


Bulk  viscosity 


1.002  -^y 


(1.16) 


Thermal  conductivity   ^     b 
/       .     ,    ,   s   A.     o 
(monatomic  molecules)   r~  -  ~ 
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-  +  1.2   +  0.757y 
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(1.17) 


where 


y  =—    X 


The  value   of  y  may  be   obtained   from   the   equation  of   state  by 


—  =   1   +  y 
RT  ^ 


(1.18) 


Enskog  suggested  as  an  empirical  modification  of  the  hard  sphere 
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theory  the  use  of  thermal  pressure  in  place  of  external  pressure. 
Thus  according  to  this  modification  y  becomes 

Equations  (1.15)  and  (1.17)  predict  that  when  (r)V)  and  (XV)  are 
plotced  as  functions  of  y,  the  curves  will  go  through  minima  at 
y  =  1.145  and  y  =  1.151  respectively,  having 

[r]V]    =  2.545r]  b^  ■    at  y  =  1.146  (1.20) 

min         "  ° 

and  [XV]    .      =  2.938X  b       at  y  =  1.151  (1.21) 

mm         o  o 

In  the  above  transport  coefficient  expressions,  b   is  the  only  unknown 

parameter.   It  may  be  obtained  from  equations  (1.20)  and  (1.21)  by 

using  experimental  values  of  [rN]    .      and  [XV]  .  . 
°        ^  '-'I   Jmin        -'mm 

The  Enskog  dense  gas  theory  gives  a  useful  approximate 
description  of  transport  coefficients  above  the  critical  temperature 
and  at  densities  less  than  the  critical  value.   The  theory  fails  at 
temperature  below  the  critical  temperature  and  at  densities  higher  than 
the  critical  density.   Under  these  conditions  the  rigid  sphere  approxi- 
mation is  in  serious  error. 

The  Enskog  theory  has  been  tested  for  thermal  conductivity  and 
viscosity  of  argon  at  various  temperatures  and  densities  (23) . 
Figure  1.3  shows  the  comparison  of  calculated  and  experimental  vis- 
cosities of  argon  at  0°C  and  75°C.   The  value  of  b  used  in  these 
calculations  was  obtained  from  equation  (1.20)  using  data  at  0  C.   The 
agreement  between  theory  and  experiment  is  moderately  good  at  0  C  for 
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densities  up  to  about  0.8  g/cc.   At  hi^er  densities  the  theory  fails 

to  predict  the  correct  density  dependence.   As  can  be  seen  in 

Figure  1.3,  the  Enskog  theory  does  not  predict  the  correct  temperature 

dependence.   Dymond  and  Alder  (25)  recently  modified  the  Enskog 

theory  by  using  temperature  dependent  rigid  sphere  diameters  obtained 

from  the  Van  der  Waals  equation  of  state.   With  this  modification,  the 

theoretical  and  experimental  values  of  viscosity  are  shown  in 

Figure  1.4.   The  predicted  temperature  dependence  is  found  to  improve 

considerably;  however  the  theory  still  fails  to  predict  the  density 

dependence  at  densities  above  0.8  g/cc. 

1 .4 Activation  and  Free  Volume  Theories 

An  activation  theory  for  liquid  transport  properties  was  first 
proposed  by  Eyring  (26).   This  type  of  theory  assumes  that  a  molecule 
spends  a  large  fraction  of  its  time  oscillating  about  an  equilibrium 
position  in  a  cell,  and  only  occasionally  does  it  leave  one  cell  to 
take  up  position  in  a  neighboring  vacant  cell.   Transport  of  mass  and 
momentum  are  assumed  to  occur  during  such  molecular  transitions. 
Later  modifications  to  Eyring' s  original  theory  have  included  the 
introduction  of  the  concept  of  the  fluctuating  free  volume.   These 
theories  are  discussed  in  detail  in  Chapter  2.   The  free  volume  theories 
apply  only  at  densities  above  approximately  twice  the  critical  value. 

Korrocks  and  McLaughlin  (27)  applied  the  activation  and  free 
volume  theory  to  the  thermal  conductivity.   They  assumed  a  face- 
centered-cubic  lattice  structure  for  the  liquid  and  that  transfer  of 
thermal  energy  down  the  temperature  gradient  was  due  to  two  causes: 
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the  actual  transit  of  a  molecule  from  one  lattice  site  to  another, 
and  that  due  to  the  collisions  of  an  oscillating  molecule  with  its 
neighbors.   For  liquids  the  contribution  due  to  the  former  cause  is 
negligible  when  compared  with  that  of  the  latter.   The  frequency  of 
oscillation  is  determined  by  the  molecular  mass  and  the  intermolecular 
force.   Their  theory  has  been  tested  for  a  number  of  simple  liquids 
and  the  agreement  between  theory  and  experiment  is  often  within  20%. 
Calculated  and  experimental  thermal  conductivities  of  a  number  of 
liquids  are  compared  in  Figure  1.5. 

1-5 The  Correspondence  Principle 

The  principle  of  corresponding  states  has  been  found  very 
useful  in  the  calculation  of  equilibrium  properties  of  dense  gases  and 
liquids  (28).   However  very  few  studies  have  been  made  of  its  applica- 
tion to  transport  properties  of  dense  gases  and  liquids.   As  the 
principle  of  corresponding  states  will  be  discussed  in  great  detail 
in  later  chapters  it  will  not  be  discussed  further  here. 
1.6     Empirical  Correlations 

A  very  good  review  of  the  purely  empirical  methods  of  estimating 
transport  coefficients  is  given  in  a  new  book  by  Reid  and  Sherwood  (29). 
These  methods  have  neither  a  theory  nor  a  model  to  describe  them.   The 
many  expressions  proposed  usually  contain  one  or  more  empirical  para- 
meters or  constants  which  are  said  to  characterize  the  structure  or 
properties  of  the  molecules;  in  most  cases  these  parameters  have  to  be 
supplied  by  the  authors  of  the  correlation.   A  typical  example  of  these 
correlations  is  Souder's  method  for  estimating  liquid  viscosity,  which 
is  one  of  the  few  empirical  methods  recommended  by  Reid  and  Sherwood. 
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Souder's  equation  is 

I. 


log(log  IOt])   =  ^  p  -  2.9  (1.22) 


where   T]  =  liquid  viscosity,  centipoise 

-3 
,0  =  liquid  density,  gem. 

I  =  viscosity  constant  calculated  from  atomic  and  structural 
constants  supplied  by  the  author 

M  =  molecular  weight 
The  results  yielded  by  such  methods  are  usually  not  very  satisfactory. 
In  spite  of  the  author's  claim  that  viscosity  predictions  are  within 
207o,  errors  are  often  greater  than  this.   Thus  the  equation  predicts  a 
viscosity  for  acetic  acid  at  40°C  that  is  36%  too  low  (29).   Similar 
empirical  expressions  are  also  available  for  liquid  thermal  conducti- 
vity and  self-diffusion  coefficient;  the  discrepancies  between  these 
equations  and  experimental  data  are  of  a  similar  order  to  that  of  the 
viscosity  correlations. 

Among  the  empirical  correlations,  those  based  on  residual- 
viscosity  concepts  have  received  much  attention  (30-32).   Figure  1.6 
is  a  plot  of  viscosity  isotherms  for  krypton  as  a  function  of  density. 
At  the  lower  densities  the  curves  are  almost  parallel,  and  this  has 
led  a  nimber  of  workers  to  correlate  the  residual  viscosity  r]  -   T]     with 
density.   The  general  form  of  this  correlation  is 

ri-ri^=  f^(p)  (1.23) 

where  f^  is  a  function  of  density  only  and  T)  is  the  dilute  gas 
viscosity.   However,  as  can  be  seen  in  Figure  1.6,  the  isotherms 
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intersect  at  higher  densities,  so  that  such  correlations  are  not  valid 
over  the  Xv^hole  density  range.   The  intersection  occurs  at  densities  of 
about  txcfice  the  critical  value.   Above  this  density  value  the  free 
volume  theory  is  found  to  work  well. 
1 .7 Summary 

Of  the  different  approaches  discussed  above,  the  time  correla- 
tion function  theory  and  the  kinetic  theory  of  liquids  are  the  most 
rigorous  and  correct  descriptions  of  transport  processes  in  dense 
fluids.   In  their  present  stage  of  development,  however,  they  offer 
little  immediate  prospect  of  yielding  methods  of  predicting  transport 
properties  for  fluids  of  interest  to  chemical  engineers. 

The  Enskog  dense  gas  theory  provides  good  results  for  dense 
gases,  but  fails  at  densities  a  little  above  the  critical  value.   For 
the  lower  densities  quite  good  results  may  also  be  obtained  from 
empirical  equations,  such  as  those  employing  the  residual  viscosity 
concept. 

At  present  a  particular  need  exists  for  satisfactory  methods  of 
predicting  transport  properties  for  fluids  at  densities  well  above  the 
critical,  that  is  in  the  normal  liquid  density  region.   Toward  this 
end  t^<ro   approaches  are  examined  in  detail  in  this  dissertation.   The 
free  volume  theory  is  first  studied,  and  examined  as  a  framework  for 
developing  predictive  corresponding  states  relations  for  simple  fluids. 
Although  such  a  model  lacks  the  desirable  rigor  present  in  the  formal 
theories,  it  has  the  considerable  advantage  of  being  solvable,  and 
suggests  ways  in  which  expressions  may  be  obtained  for  mixtures. 
Chapter  2  briefly  reviews  previous  work  on  free  volume  theories,  and 
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Chapters  3  and  4  contain  the  new  contributions  to  the  theory. 

The  second  approach  consists  of  the  development  of  corresponding 
states  relationships  directly  from  the  time  correlation  function  theory. 
The  correspondence  principle  is  particularly  valuable  when  based  on 
the  rigorous  statistical  mechanical  approach,  but  no  serious  attempt  to 
apply  che  principle  to  dense  fluid  transport  properties  seems  to  have 
been  made  previously.   This  approach  is  examined  in  detail  in  Chapters 
5  and  5 . 


CHAPTER  2 

PREVIOUS  WORK  ON  ACTIVATION  AND 
FREE  VOLUME  THEORIES 


Because  of  the  present  difficulties  associated  with  developing 
a  rigorous  kinetic  theory  of  dense  fluids  analogous  to  that  available 
for  dilute  gases,  considerable  attention  has  been  given  to  simplified 
models  of  the  liquid  state.   Assumptions  are  introduced  in  the  initial 
stages  of  the  development,  an  attempt  being  made  to  obtain  a  model 
which  incorporates  the  essential  features  of  real  liquids,  but  which 
is  still  solvable.   Such  attempts  are  exemplified  by  the  activation  and 
free  volume  theories  for  liquid  transport  properties. 

Eyring  (26)  may  be  considered  a  pioneer  in  the  development  of 
activation  and  free  volume  theories  for  liquid  transport  properties. 
In  Eyring' s  original  approach  the  equilibrium  positions  of  molecules 
in  a  liquid  were  considered  to  be  on  a  regular  lattice  similar  to  the 
crystal  lattice  in  a  solid.   The  transport  of  mass  and  momentiim  occurs 
as  a  result  of  a  molecule  jumping  from  one  equilibriiom  lattice  site 
onto  a  vacant  neighboring  site.   In  order  to  make  such  a  jxjmp  a 
molecule  is  required  to  have  the  necessary  activation  energy  to  over- 
come the  minimum  energy  barrier  separating  two  adjacent  sites,  and  at 
the  same  time  a  vacant  site  must  be  available.   Using  this  model  Eyring 
obtained  the  following  expressions  for  liquid  viscosity  and  self- 
diffusivity  by  making  use  of  the  absolute  reaction  rate  theory: 

1    ,^   ,^xl/2  3  %^^^ 


T]   =  -fY3(271inkT)"'  q^e  "  (2.1) 

pv 
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2  -e^/kT 

^  =  1/2   1/3  ^  ^2-^> 

(27nnkT)^^  q^  '^ 

More  recently  Weymann  (35,36)  used  a  statistical  mechanical 
approach,  and  arrived  at  expressions  very  similar  to  those  of  Eyring. 
In  his  derivation  of  the  equations  for  transport  properties,  Weymann 
gave  a  clearer  picture  of  the  physical  model  used  and  of  the  signifi- 
cance of  the  activation  energy.   Moreover  his  approach  suggests  a 
possible  extension  to  liquid  mixtures.   Weymann  also  derived  an  expres- 
sion in  which  the  probability  of  hole  formation  is  related  to  the 
volume  and  energy  needed  to  form  a  hole.   The  equations  obtained  by 
Weymann  for  viscosity  and  self-diffusion  coefficient  are: 


^  =  4^f  "in)  (2TltnkT)^/2e'°  (2.3) 


2X,  1  .  .,  „  V.  1/2   -e^/kT 
D  =— ll-   2kT  \     ^  o  (2.4) 

3   \   Tim  / 


where  Xy,   is  the  probability  of  hole  formation. 

Eyring' s  activation  theory  was  modified  by  McLaughlin  (37) 
by  introducing  the  concept  of  the  probability  of  hole  formation,  as 
was  done  by  Weymann.   In  his  treatment,  he  assumed  a  face-centered- 
cubic  lattice  for  the  liquid.  Tae   most  important  modification  made  by 
McLaughlin  was  the  attempt  to  relate  the  activation  energy  and  energy 
of  hole  formation  to  the  intermolecular  forces  by  means  of  the 
Lennard- Jones  and  Devonshire  cell  theory  of  liquids.   A  similar 
expression  for  liquid  viscosity  was  also  obtained  by  Majumdar  (38) 
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using  the  tunnel  model  of  Barker  (39)  in  place  of  the  Lennard- Jones 

and  Devonshire  cell  model. 

The  Eyring  and  Weymann  theories  yield  good  results  for  the 

viscosity  of  Arrhenius  liquids  (ones  whose  viscosity  varies  with 

temperature  according  to  an  equation  of  the  type  Inr)  =  C  +  B  /T)  , 

a    a 

but  cannot  account  for  non- Arrhenius  liquid  behavior.   In  order  to 
overcome  this  deficiency  Doolittle  (40-42)  proposed  empirical  expres- 
sions which  related  liquid  transport  properties  to  the  free  volume, 
defined  by 

v^  =  V  -  V  (2.5) 

f        o 

where  v  is  the  hard-core  volume  of  the  molecule.   His  free  volume 
o 

equations  were  placed  on  more  solid  theoretical  grounds  by  Cohen  and 

Turnbuli  (43,44)  who  derived  them  by  a  statistical  mechanical  method. 

They  assumed  that  a  molecule  moves  about  in  a  cell  in  a  gas- like 

manner,  while  the  free  volume  available  to  each  molecule  fluctuates 

with  time.   A  molecule  is  able  to  jump  into  a  neighboring  cell  only  if 

a  free  volume  of  a  size  greater  than  some  minimum  value  v   is  available. 

°  o 

They  obtained  the  following  expressions  for  the  viscosity  and  self- 
diffusion  coefficient: 


Tne  term  "free  volume"  in  these  theories  has  a  different 
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meaning  from  that  implied  in  thermodynamic  free  volxome  theories.   The 
theories  of  Doolittle  and  Cohen  and  Turnbull  assumed  the  molecules  to 
behave  as  hard  spheres;  the  free  volume  referred  to  here  is  the  space 
in  the  fluid  unoccupied  by  the  hard  spheres  themselves.   The  average 
free  volume  per  molecule  is  defined  as  v  -  v  ,  where  v  =  V/N.   The 
free  volume  theory  of  these  authors  describes  the  viscosity  behavior 
at  atmospheric  pressure,  but  fails  to  predict  the  temperature  depen- 
dence of  viscosity  at  constant  volume. 

Recently,  Macedo  and  Litovitz  (45)  proposed  expressions  for 
viscosity  and  self-diffusion  coefficient  in  which  the  hole  probability 
in  Weymann's  equations  was  replaced  by  the  free  volume  expression 
obtained  by  Cohen  and  Turnbull.   Molecular  transport  is  assumed  to 
occur  if  a  molecule  has  sufficient  energy,  e   to  overcome  intermolecular 
forces  with  its  neighbors,  and  at  the  same  time  it  has  a  free  volume 
greater  than  the  minimum  free  volume  v  needed  for  a  jump  to  occur. 
The  equations  proposed  for  viscosity  and  self-diffusion  coefficient 


are 


7V 
T]  =   Aexp  f  -—V  )  exp  f^  )  (2.8) 


V  -  V    ;     '^  VRT 

o 


7V    N.       /  E  \ 

D  =  Bexp  (    -   —V  ]  exp  (-   ^  )         (2.9) 


V  -  V   ;   "     \      KT 


J 


By  treatirg  V  and  E   as  empirical  constants  Macedo  and  Litovitz 
have  shown  that  equation  (2.8)  describes  the  viscosity  behavior  of  a 
number  of  liquids  over  a  range  of  temperature.   Chung  (46)  has  pre- 
sented an  elegant  statistical  mechanical  derivation  of  the  equation  of 
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Macedo  and  Litovitz.   His  derivation  is  presented  in  Appendix  1. 

Macedo  and  Litovitz  assume  the  pre-exponential  factor  A  in 

equation  (2.8)  to  be  proportional  to  temperature  T  although  most  other 

1/2 
vjorkers  predict  a  temperature  dependence  of  T    .   Both  types  of 

temperature  dependence  of  the  pre-exponential  factor  have  been  tested 

in  this  work,  and  experimental  results  seem  to  give  better  agreement 

1/2 
with  theory  xvhen  a  temperature  dependence  of  T    is  used.   Thus 

throughout  this  work  the  pre-exponential  factors  A  and  B  are  assumed 

1/2 
to  be  proportional  to  T  '  .   Thus  equations  (2.8)  and  (2.9)  may  be 

rewritten  as 


(2.10) 


where 


and 


V 

=  A  exp  / 

^       ^o       N             / 

,1/2 

.'-\J^^'K 

RT    I 

D 

=  B^exp  ( 

^    \    \ 

/     E 
1         ^ 

,1/2 

V  -  V     I-' 

^                   o 

y     RT 

A      = 
o 

A 
,1/2 

B      = 
o 

B 
.1/2 

(2.11) 


T 
Macedo  and  Litovitz  original  equations  contain  a  constant 

factor  7  which  was  introduced  to  account  for  the  overlapping  of  free 

volumes.   In  the  early  part  of  this  work,  the  constant  7  was  evaluated 

for  a  number  of  molecules  and  x^7as  found  to  be  close  to  unity.   This 

factor  will  therefore  not  be  included  in  the  equations  of  viscosity  and 

dif fusivity . 


CHAPTER  3 
IMPROVED  FREE  VOLIME  THEORY 

The  viscosity  equation  (2.8)  has  been  tested  by  Macedo  and 
Litovitz  (45)  for  a  variety  of  liquids  and  by  Kaelble  (47)  for  poly- 
meric substances.   They  treated  the  pre-exponential  factor  A  ,  the 

activation  energy  E   and  the  minimum  free  volume  V  as  adiustable 
"■^   V  o      -^ 

parameters  characteristic  of  the  substance  considered.   Equation 

(2.8)  was  found  to  describe  the  viscosity  behavior  of  the  substances 

tested  by  these  workers  moderately  well.   On  the  other  hand,  Naghizadeh 

and  Rice  (18)  tested  the  theory  for  the  sclf-dif fusivity  of  simple 

fluids  (such  as  the  inert  gases)  and  found  that  the  agreement  between 

theory  and  experiment  was  poor,  especially  in  the  high  density  region. 

If  V   is  a  constant,  as  assumed  by  Macedo  and  Litovitz, 
o  ^ 

differentiation  of  equation  (2.10)  with  respect  to  1/T  at  constant 
volume  (i.e.  constant  density)  yields 

5ln(Ti/T^^^  1   =^ 

a(l/T)    J^    R  (3.1) 

E  represents  the  minimum  energy  required  by  a  molecule  to  overcome 
intermolecular  forces  in  making  a  jump  and  was  also  assumed  constant 

by  Macedo  and  Litovitz.   Therefore  equation  (3.1)  predicts  that  a 

1  /  ? 
plot  of  ln('r)/T  '  )  versus  reciprocal  temperature  at  constant  volume 

should  give  a  series  of  parallel  straight  lines.   For  a  variety  of 

nonpolar  and  slightly  polar  liquids  over  a  moderate  density  and 

temperature  range  such  a  plot  produces  straight  lines,  but  the  slopes 


27 


28 


of  the  lines  vary  with  volume.   Moreover,  when  the  results  are  plotted 
over  a  wide  temperature  range  nonlinearity  becomes  apparent.   Such 
nonlinearity  is  particularly  marked  for  fluids  composed  of  simple 
molecules.   This  is  illustrated  in  Figures  3.1  and  3.2  for  argon  and 
nitrogen,  for  which  data  are  available  at  constant  volume  over  wide 
ranges  of  temperature  and  density. 

From  the  above  discussion  it  is  apparent  that  the  equations  as 
used  by  Macedo  and  Litovitz  do  not  correctly  predict  qualitatively 
the  effect  of  temperature  and  volume  on  the  viscosity,  especially  for 
simple  fluids.   In  addition,  extensive  experimental  viscosity  data  are 
needed  for  each  fluid  in  order  to  fit  the  adjustable  parameters  A  , 
E   and  V  :  when  used  in  this  way  their  equation  is  no  more  than  an 
empirical  correlation. 

In  this  chapter  the  physical  significance  of  the  parameters 
E   and  V  is  examined  in  the  light  of  the  theory,  and  their  dependence 
on  molecular  type,  temperature  and  volume  is  discussed.   Interpretation 
of  these  quantities  on  the  molecular  level  suggests  corresponding 
states  relationships  which  may  be  used  to  predict  the  parameters.   The 
improved  theory  also  explains  why  the  theory  of  Macedo  and  Litovitz 
fails  for  simple  fluids. 
3.1     Temperature  Dependence  of  V 

The  parameter  V  of  equations  (2.10)  and  (2.11)  represents  the 
minimum  free  volume  that  must  be  available  before  a  jump  may  occur.   If 
the  molecules  may  be  treated  as  rigid  spheres,  as  in  the  smoothed 
potential  cell  model  of  Prigogine  (50),  the  free  volvime  is  independent 
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Figure   3.1.      Viscosity   of  Argon. 
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of  the  teniDerature  at  constant  volume,  and  V  should  therefore  be  a 

o 

constant  for  a  particular  molecule,  independent  of  temperature  and 

density.   Macedo  and  Litovitz  (45)  assumed  this  to  be  the  case  for  all 

molecules  treated  by  them. 

For  more  realistic  potential  models  one  would  expect  V   to 

o 

decrease  somewhat  with  increasing  temperature,  since  as  temperature 
rises  the  average  kinetic  energy  of  the  molecules  increases,  and 
molecules  are  thus  able  to  approach  each  other  more  closely.   According 
to  such  a  viewpoint  V   should  not  be  affected  by  the  density  at  constant 
temperature.   The  extent  to  which  this  parameter  depends  on  temperature 
will  be  determined  largely  by  the  repulsive  portion  of  the  intermolecu- 
lar  potential  energy  curve.   Since  this  part  of  the  curve  rises  less 
steeply  for  simple  molecules  than  for  more  complex  polyatomic  molecules, 
one  would  expect  the  effect  of  varying  V   to  be  most  evident  for  the 
simpler  molecules. 

To  obtain  a  general  expression  describing  the  temperature 
dependence  of  V  ,  a  [6,n]  pair  potential  energy  function  is  assumed. 
Tnus 


a  \  fa 


-I 


(3.2) 


u  =  £K 


w^s^e  /       \   X   \  6/(n-6) 

K  =  ^   n   \  /n^ 
1 


The  parameter  n  indicates  the  steepness  of  the  repulsive  part  of  the 
curve;  for  small  molecules  such  as  the  inert  gases  n  is  close  to  12, 
whereas  for  more  complex  molecules  the  best  values  of  n  is  28  or 
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higher. 

For  a  group  of  roughly  spherical  molecules  V  may  be  assumed 

to  be  proportional  to  the  cube  of  the  molecular  "diameter."  An 

estimate  of  the  variation  of  the  molecular  diameter  (and  hence  V  )  with 

^  o' 

temperature  may  be  obtained  by  equating  the  average  kinetic  energy  of 
a  t^v^o-particle  system  to  the  potential  energy  of  the  system  at  the 
distance  of  closest  approach,  d: 


3kT  =  £K, 


(3.3) 


or,  in  reduced  form 


3T*  =  K, 


n/3 


2  -1 


(3.4) 


where 


T""  = 


kT 


V 


Equation  (3.4)  may  be  solved  to  obtain  v  '  as  a  function  of  T'  for 

a 

various  n  values.   The  temperature  dependence  of  v'!  for  several  n 
values  is  shown  in  Figure  3.3.   Assuming  that  V  is  proportional  to 
d  ,  one  may  write 


V  =  ^v 
o    ^   a 


(3.5) 


where  ^  is  constant  for  a  particular  molecule,  V  at  any  temperature 
can  be  calculated  from 


V 


V    =  V    •  -rr S: 

o     o    V"  (T"  =  1) 

a 


(3.6) 
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where  V  is  the  value  of  V  at  T'  =  1.   Thus  a  knowledge  of  V  , 
o  o  o 

together  with  n  and  e,  suffices  to  calculate  V  at  any  temperature. 
Some  support  for  the  above  procedure  is  provided  by  values  of  the  hard- 
sphere  diameter  calculated  at  various  temperatures  from  equation  of 
state  data  for  inert  gases  by  Dymond  and  Alder  (25) .   These  values  are 
included  in  Figure  3.3,  and  agree  well  with  the  curve  for  n  =  12. 

With  n  =  coj  no  temperature  dependence  should  be  observed  for 
V  .   Under  these  conditions,  from  equation  (2.10) 

'   5ln(VT^^Sl   _  ^v  ,^   ,. 

S(l/T)   J^  -"i  ^^-^^ 

1/2  - 1 

•When  ln(rj/T  '  )  is  plotted  against  T   at  constant  volume  for  such 

nonpolar  molecules  as  decahydronaphthalene,  benzene,  etc.,  approximately 
linear  behavior  is  observed;  for  these  molecules  n  is  large  and  the 
variation  of  V  with  temperature  is  small.   For  simple  molecules,  such 
as  the  inert  gases,  for  which  n  %  12,  the  temperature  dependence  is 
more  marked,  and  the  model  predicts  noticeable  nonlinearity  on  such  a 
plot.   This  is  as  observed  experimentally,  as  shown  in  Figure  4.1  in 
the  next  chapter. 

In  addition  to  the  above  comments,  the  parameter  V  should  be 
related  in  some  way  to  molecular  size.   If  V  may  be  made  dimension- 
less  with  some  suitable  reducing  parameter,  the  resulting  reduced 
minimum  free  volume  may  be  quite  constant  for  a  series  of  similar 
molecules.   Since  free  volume  theories  apply  best  at  high  densities, 
the  molal  volume  at  the  melting  point  seems  a  more  appropriate  reducing 
parameter  than  the  critical  volume.   Thus 
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1-    '' 

,,i"      o 

^o   =F  (3.8) 

m 


3.2     Volume  Dependence  of  E 

Macedo  and  Litovitz  (45)  treated  the  minimum  activation  energy 
E   as  an  empirical  constant,  independent  of  density  and  temperature. 
However,  as  pointed  out  by  Brummer  (51)  ,  E  may  be  expected  to  vary 
with  the  average  intermolecular  distance,  and  thus  with  density.   The 
value  of  E  becomes  larger  as  the  molal  volume  decreases  because  of 
the  increase  in  repulsive  intermolecular  force  between  molecules  at 
close  separations. 

From  equations  (2.10),  (3.4)  and  (3.6)  one  may  obtain 


E    ,. 
V  _  Lim 

R   T_».0 


^inlnZl^l  (3.9) 

a(l/T) 


J  V 


so  that  E  may  be  calculated  from  experimental  viscosity  values.   For 

1/2 
more  complex  molecules  where  n  is  large,  a  plot  of  ln('n/T   )  versus 

1/T  is  found  to  be  approximately  linear,  and  the  requirement  T  — ♦■  0  in 
equation  (3.9)  is  less  stringent.  Figure  3.4  shows  the  volume  depen- 
dence of  E   for  several  fluids  as  calculated  from  equation  (3.9). 

The  activation  energy  arises  from  the  motion  of  the  jumping 
molecule  from  its  initial  equilibrium  position  through  a  region  of 
higher  potential  energy  to  its  final  position.   A  model  of  the  situation 
is  shown  in  Figure  3.5  in  which  a  molecule  jumps  from  A  to  C ,  and 
passes  through  a  region  B  in  which  it  has  to  squeeze  through  a  ring  of 
2  molecules.   The  activation  energy  may  be  written 
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Figure  3,5.   Model  of  a  Molecular  Jump. 
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Define 


e^  =  z  [uCr^)  -  u(r^) 


=  K  Z£ 


r  /    \  n 


c,  =  v/r: 


a_ 


K  Z£ 


0_ 
r. 


(3.10) 


c_  =  c 


(r2/ri) 


C/.  = 


kT 
m 


Then 


e   =  K.zkT  c, 
V    1   m  4 


I       3  n 


i  c^v 


n/3 


3  m 


(3.11) 


The  parameters  c^,  z^,    c,,  c^  and  n,  the  repulsive  ej/ponent  should  be 
approximately  the  same  for  a  group  of  similar  molecules  and  thus 
equation  (3.11)  suggests  a  corresponding  states  relation  of  the  form 


E"  =  E"(V'') 
V     V    "^ 


(3.12) 


where 
and 


V    RT 


m 


r^,  V_ 

V'"  =  V 


m 


Figure  3.6  shows  the  correlations  of  reduced  activation  energy  as  a 
function  of  reduced  molal  volume  for  several  nonpolar  molecules.   Good 
agreement  is  obtained  at  high  reduced  volumes,  but  some  scatter  is 
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observed  at  lower  values  of  v",  where  the  curve  rises  steeply.   This 
may  be  attributed  to  differences  in  n,  the  repulsive  exponent,  for  the 
molecules  within  the  group.   Thus  t\-Jo   curves  A  and  B  are  plotted  in 
Figure  3.6.   Curve  A  is  for  simple  molecules  obeying  a  Lennard- Jones 
[6,12]  potential.   Curve  B  shows  an  approximate  relationship  for  more 
complex  molecules  which  obey  a  potential  law  in  which  n  is  larger 
than  12.   Sources  of  viscosity  data  from  which  values  of  E   are  cal- 
culated are  given  in  Figure  3.6. 

Taking  into  consideration  the  temperature  dependence  of  V^ 
and  the  volume  dependence  of  E  ,  the  modified  viscosity  and  self- 
diffusion  coefficient  equations  become: 


.   V  (T)   X      .  E  (V)  X 
^     -'    -/^— ^ Uxp(^^^  (3.13) 


TT^-  =  A  exp  f 
T^^^     °    ^-  V  -  V  (T) 


and 


/    V  (T)   \     /   E  (V) 
°    =  B  exp  (-   T7-Vt^  )  exp  r  -  ^—  ]  (3.14) 


^1/2    o""^  \^  V  -  V^(T)  J  '-^'^  \  RT 


CHAPTER  4 
TEST  OF  IMPROVED  FREE  VOLUME  THEORY 

In  order  to  perform  a  rigorous  test  of  the  theory,  viscosity 
and  self-diffusion  coefficient  data  over  a  wide  range  of  temperature 
and  density  conditions  are  most  desirable.   In  this  work  only  fluids 
composed  of  simple  molecules  which  approximately  obey  the  Lennard- Jones 
(6,12)  potential  law  will  be  considered.   These  fluids  include  the 
inert  gases  and  most  diatomic  fluids- and  perhaps  methane,  fluoro- 
methane  and  te trade uterome thane. 

The  only  monatomic  and  diatomic  fluids  for  which  viscosity 
data  over  a  wide  range  of  temperature  and  density  conditions  are  avail- 
able seem  to  be  argon  and  nitrogen.  Zhadanova  (48,49)  has  reported 
viscosity  measurements  for  these  fluids  for  experimental  conditions 
corresponding  to  pressures  up  to  several  thousand  atmospheres.   The 
accuracy  of  the  experimental  argon  data  appears  to  be  of  the  order  67o 
over  the  entire  range  of  conditions.   However,  the  nitrogen  data  are 
in  poorer  agreement  with  measurements  of  other  workers,  and  the  accuracy 
appears  to  be  of  the  order  12-157o. 

The  only  high  pressure  self-dif fusivity  measurements  for  simple 
liquids  seem  to  be  those  of  Naghizadeh  and  Rice  (18) ,  who  reported 
measurements  for  argon,  krypton,  xenon  and  methane  over  a  reasonable 
range  of  temperature  and  at  pressures  up  to  a  more  than  one  hundred 
atmospheres.   Naghizadeh  and  Rice  claimed  that  the  uncertainty  of 
their  experimental  data  was  less  than  5%.      However  the  disagreement 
between  their  argon  data  and  those  of  Corbett  and  Wang  (56)  is  about 
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12?o,  while  their  xenon  data  differ  from  those  of  Yen  and  Norberg 
(57)  by  as  much  as  86%  and  their  reported  methane  data  are  about  11% 

higher  than  those  of  Gaven,  Waugh  and  Stockmayer  (58). 

4.1 Corresponding  States  Relationships 

The  following  procedures  were  employed  to  test  the  theory  for 

viscosity  and  to  establish  the  corresponding  states  relationships  in 

equations  (3.8)  and  (3.12). 

(1)  The  experimental  data  were  plotted  as  ln(r]/T^''^)    against 

T   at  constant  volim:ie.   Values  of  E^  were  estimated  at  various  volumes 
from  the  lines  extrapolated  to  low  temperature,  since  from  equation 
(3.9) 

!v  ^  Lim  r  ^IM-n/T'^^)    1 
R   T->0  [_    8(1/T)   J  ^^■■^^ 

(2)  Values  of  the  parameters  A  and  V  were  obtained  from 

the  best  fit  to  the  viscosity  data  using  the  computer,  equations  (3.4) 

and  (3.5)  being  usee  to  obtain  the  temperature  deoendence  of  V 

'  o 

Values  of  E^  for  argon  and  nitrogen  estimated  by  procedure  (1) 
are  included  in  Figure  3.6.   Tneoreticai  and  experimental  liquid  argon 
viscosity  values  are  shown  in  Figure  4.1.   Agreement  between  theory  and 
experiment  is  found  to  be  within  5%.   Similar  agreement  between  theory 
and  experiment  is  found  for  nitrogen,  with  V"""  =  17.0  cc/mole  and 
A^  =  2.53x10   poise  (°K)    ~,  and  is  shown  in  Appendix  2.   The  theory 
provides  a  satisfactory  fit  to  the  data  for  densities  above  twice  the 
critical  value.   Thus  as  can  be  seen  in  Figure  4.1  the  theory  fails 
for  molal  volumes  of  39.15  cc/mole  and  above  for  argon.   Similar 
behavior  is  also  observed  for  the  nitrogen  viscosity  (see  Figure  A2.1 
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in  Appendix  2),  where  the  theory  fails  for  molal  volumes  of  56.0  cc/mole 
and  above.   The  density  range  in  which  the  theory  applies  is  sharply 
defined  and  is  similar  to  that  in  which  similar  models  for  thermo- 
dynamic properties,  such  as  the  Prigogine  smoothed  potential  model 
(50),  apply.   The  theory  appears  to  work  well  for  temperatures  below 
about  1.5  T  ,  although  the  temperature  range  in  which  the  theory 
applies  is  less  clearly  defined. 

The  values  of  the  parameters  E  and  V   found  for  argon  and 

^         .  V      o  ° 

nitrogen  may  be  used  to  form  the  basis  of  corresponding  states  cor- 
relations of  these  parameters  as  described  in  Chapter  3,  from  which 
viscosity  and  self-dif fusivity  values  of  other  molecules  may  be  pre- 
dicted.  The  average  value  of  the  reduced  minimum  free  volume  for 
argon  and  nitrogen  is  0.554,  so  that 

vj"  =^  =  0.554  (4.2) 

m 


L    9 


Test  of  Proposed  Correlations 


To  use  equations  (2.10)  and  (2.11)  to  estimate  the  viscosity 

and  self-diffusion  coefficient  for  simple  molecules,  values  of  V  were 

o 

calculated  from  equations  (3.4),  (3.6)  and  (4.2),  assuming  n  =  12  and 
taking  values  of  e/k  determined  from  gas  viscosity  data  (20) .  Values 
of  E  were  found  using  curve  A  of  Figure  3.6.  The  best  values  of  A 
were  obtained  by  substituting  experimental  viscosity  data  in  equation 
(2.10)  for  each  molecule.  In  a  similar  manner,  by  fitting  experimental 
self-diffusion  coefficient  data  of  each  molecule  into  equation  (2.11) 
the  best  values  of  B  were  obtained.   Other  than  the  viscosity  data 
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for  argon  and  nitrogen  of  Zhadanova  (48,  49)  and  self-diffusion 
coefficient  data  of  Naghizadeh  and  Rice  (18),  very  few  high  pressure 
viscosity  and  self-dif fusivity  data  are  available.   Thus  most  of  the 
viscosity  and  self-dif fusivity  data  used  are  for  saturated  liquids. 
Table  4.1  shows  the  viscosity  data  sources  and  ranges  of  temperature 
and  density  for  each  of  the  molecules  studied.   The  different  para- 
meters which  are  required  for  the  estimation  of  viscosity  with 
equation  (2.10)  and  the  average  percent  deviation  between  theory  and 
experiment  are  shovm  in  Table  4.2.   Tables  4.3  and  4.4  provide 
similar  information  for  the  self-diffusion  coefficient  of  the  molecules 
studied. 

The  average  percent  deviation  bettvreen  theory  and  experiment 
is  about  67o  or  less  for  viscosity.   For  self-diffusion  coefficient  the 
discrepancy  between  experiment  and  theory  is  found  to  be  much  larger, 
being  about  157o  in  most  cases.   The  discrepancies  between  theory  and 
experiment  for  both  viscosity  and  self-diffusion  coefficient  seem  to 
be  of  the  same  order  as  the  accuracy  of  che  experimental  data  for  the 
liquids  studied.   The  viscosity  data  of  Zhdanova  for  nitrogen  at 
high  pressures  seem  to  be  in  poor  agreement  with  those  of  other  workers 
(66),  and  saturated  liquid  viscosity  data  (66)  were  used  in  obtaining 
the  Aq  value  given  in  Table  4.2  for  this  fluid. 

Experimental  and  predicted  viscosities  for  xenon  and  methane 
are  compared  in  Figures  4.2  and  4.3  (similar  plots  for  Ne,  Kr  and  N 
are  shovjn  in  Appendix  2)  .   Tne  predicted  values  become  less  reliable 
as  density  falls  to  values  approaching  Ip     (for  methane  this  corresponds 
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Free  Volume  Theory 

o    Experiment  (64,65) 


•7.5 


-8.0  j- 
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//  R   3^-1 
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Figure  4.2.   Test  of  Free  Volume  Theory  for 
Liquid  Xenon  Viscosity. 
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to  V  =  49.5  cc/mole) .   Figure  4.4  compares  theory  and  experiment  for 

the  self-diffusion  coefficient  of  saturated  liquid  methane. 

The  expressions  for  the  pre-exponential  factors  of  the 

viscosity  and  self-diffusion  coefficient  equations  of  Eyring  (26) 

2/3  ,,,1/2 


and  Weymann  (35)  suggest  that  the  groups  (A^V    /M   )  and 

(B  M   /V    )  may  be  approximately  the  same  for  different  molecules. 

2/3   1/2 

Values  shown  in  Table  4.2  indicate  that  the  quantity  (A  V    /M   ) 

is  approximately  constant  for  inert  gas  liquids,  but  for  the  other 

molecules  a  range  of  values  is  found.   Similar  behavior  is  observed 

1/2   1/3 
for  the  group  (B  M   /V^        ) .   Apparently  the  expressions  given  for 

the  pre-exponential  factors  by  Eyring  and  Weymann  are  not  generally 

valid,  and  cannot  form  the  basis  of  successful  corresponding  states 

relationships  for  A^  and  B  .   In  general  the  value  A   (or  B  )  can  be 

determined  from  a  single  experimental  value  of  viscosity  (or  self- 

diffusivity)  for  a  given  fluid.   This  can  then  be  used  together  with 

the  proposed  correlations  for  V  and  E   to  predict  viscosity  (or  self- 

diffusivity)  values  at  other  temperatures  and  densities. 

4.3 Summary 

By  allowing  for  the  variation  of  V  with  temperature  it  is 

possible  to  account  quantitatively  for  the  nonlinearity  of  constant 

volume  plots  of  ln(ryT-'"'^)  vs.  T"-*"  and  ln(D/T-'-''^)  vs.  T"  "'■  for  simple 

liquids.   The  physical  model  predicts  a  density  dependence  of  E  ,  and 

thus  explains  the  variation  with  density  of  the  slope  of  constant 

volume  plots  of  ln(T]/T  '  )  vs.  t'   and  ln(D/T-'-/^)  vs.  t'''".   With  these 

improvements  the  free  volume  theory  is  able  to  accurately  describe 

the  viscosity  and  self-dif fusivity  behavior  of  liquids  composed  of 


53 


•a 

•e-l 

3 


u 

3 

•w 

« 

M 

iJ 

O 

IM 

>S 

l-J 

• 

O 

>^ 

01 

•u 

J= 

•H 

H 

> 

•H 

0) 

W 

B 

^ 

3 

liJ 

1— 1 

U-l 

O 

•H 

> 

o 

<u 

14-1 

(U 

I-H 

u 

4) 

fu 

M 

iw 

(U 

o 

C 

n 

u 

x; 

CO 

ij 

Oi 

e) 

H 

S 

<■ 

vf 

<u 

u 

3 

00 

I-  Z      ^ 


54 


simple,  nonpolar  molecules  that  may  be  expected  to  approximately  obey 
a  Lennard- Jones  [6,12]  potential  law.   The  theory  works  well  at 
densities  p  >  2^  ,  where  other  correlations  fail. 

The  fact  that  both  V  and  E  values  vary  in  a  corresponding 
states  way  indicates  that  the  free  volume  theory  provides  a  reasonably 
correct  picture  of  viscosity  and  self-dif fusivity  behavior  for  these 
liquids,  and  shows  that  the  equations'  success  does  not  arise  only 
from  the  availability  of  three  adjustable  parameters.   The  simple 
expressions  given  for  A  and  B   by  Eyring  and  others,  however,  do 
not  seem  to  be  generally  applicable. 


CHAPTER  5 

PRINCIPLE  OF  CORRESPONDING  STATES  FOR 
MONATOMIC  FLUIDS 


The  free  volume  equations  discussed  in  Chapters  2-4  accurately 
predict  viscosities  and  self-diffusion  coefficients  for  nonpolar 
fluids  in  the  density  range  p  >  2p  .   However,  the  theory  cannot  be 
readily  extended  to  include  thermal  conductivity,  and  is  limited  to 
the  above  density  range.   In  this  chapter  and  the  next  a  general  and 
direct  application  of  the  corresponding  states  principle  is  examined. 
The  treatment  for  monatomic  molecules  does  not  assume  any  simple 
model,  and  should  be  applicable  to  all  transport  properties  over  the 
entire  range  of  densities  and  temperatures. 
5 . 1 Molecular  Basis  of  the  Correspondence  Principle 

If  correctly  applied,  the  principle  of  corresponding  states 
provides  an  accurate  and  very  useful  method  for  predicting  physical 
properties  of  pure  substances,  under  conditions  for  which  no  data  are 
available,  from  the  measured  properties  of  one  or  more  substances. 
The  corresponding  states  principle  has  been  widely  applied  to  the 
thermodynamic  properties  of  pure  substances;  thus  the  generalized 
charts  of  Hougen,  Watson  and  Ragatz  (28)  and  of  Hirschfelder ,  Curtiss 
and  Bird  (20)  have  proven  of  great  practical  use  in  engineering 
calculations.   The  principle  of  corresponding  states  was  first 
suggested  by  Van  der  Waals  in  1873  as  a  result  of  his  equation  of 
state.   It  was  not  until  1939,  however,  that  a  rigorous  and  more 
general  theoretical  derivation  of  the  principle  was  attempted  by 
Pitzer  (74)  for  thermodynamic  properties  using  statistical  mechanics. 
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Pitzer  presented  his  derivation  based  on  the  following  assimptions; 

Assumption  1.   Quantiim  effects  are  negligible. 

Assimiption  2.   The  intramolecular  degrees  of  freedom 

(rotational  and  vibrational)  are  indepen- 
dent of  density. 


Assumption  3, 


The  intermolecular  potentials  are  pair- 
wise  additive. 


Assumption  4.   The  potential  energy  for  a  pair  of  molecules 
has  the  form 


u  =  eg 


(5.1) 


where   r  =  the  intermolecular  distance 
£  =  characteristic  energy 
a  =   characteristic  distance 
0  =  a  universal  function 

Using  Assumptions  3  and  4  the  configuration  energy  may  be 
written  as 


N 


U  = 


a 


(5.2) 


i<j 


The  configurational  partition  function  may  then  be  put  in 
reduced  form  (50) 


-■h 


exp 


0 


ii^ 


^  J  i 


dr^  ....  dr^^ 
—1       — N 


(5.3) 


z  =|f  z''(T^v^N) 


(5.4) 
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where   N  =  number  of  molecules 

•A-   kT 

T  =  ~~  =  reduced  temperature 

V 
V"  =  ~~T  =  reduced  molal  volume 

a 
Since  all  equilibrium  conf igurational  properties  of  interest 

may  be  evaluated  from  the  configurational  partition  function  Z  and 
its  temperature  and  volume  derivatives,  equation  (5.4)  represents  a 
general  proof  of  the  principle  of  corresponding  states  for  thermo- 
dynamic properties  of  pure  substances,  provided  all  the  above 
assumptions  are  observed.   Thus  the  reduced  equation  of  state  is 


e         e  V   ov  , 


=  P"(T",V")  (5.5) 

where  it  is  noted  that  intensive  properties  do  not  depend  upon  N. 
The  principle  of  corresponding  states  for  thermodynamic  properties  may 
also  be  derived  by  the  method  of  dimensional  analysis  (75) .   However 
in  this  method  the  conditions  under  which  a  molecule  may  obey  the 
corresponding  states  principle  are  not  explicit. 

One  way  of  deriving  the  principle  for  transport  properties 
is  by  dimensional  analysis,  with  the  reduction  performed  using  either 
the  critical  constants  or  appropriate  combinations  of  molecular 
parameters.  A  rigorous  statistical  mechanical  derivation  for  trans- 
port properties  similar  to  that  mentioned  above  for  thermodynamic 
properties,  is  possible, as  has  been  shown  by  Helfand  and  Rice  (76). 
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The  transport  coefficients  may  be  expressed  in  terms  of  equilibrium 
time  correlation  functions.   The  equations  for  the  shear  and  bulk 
viscosity,  thermal  conductivity  and  self-diffusion  coefficients  are 
given  in  equations  (1.1)  to  (1.4)  of  Chapter  1.   Using  the  four 
assumptions  made  in  deriving  the  principle  for  thermodynamic  pro- 
perties, it  is  possible  to  put  these  equations  in  a  reduced  form. 
The  derivation  of  the  principle  for  the  self-diffusion  coefficient  is 
given  below  as  an  example  of  the  procedure. 

The  self-diffusion  coefficient  expressed  in  terms  of  the  time 
integrals  of  appropriate  autocorrelation  function  has  the  form 


D  =   /  <V^(0)V^(t)>  dt  (5.6) 

J    X     X 

o 


Based  on  Assumptions  1  and  2  given  above  the  time  correlation  function 
<V^(0)V^(t)>  may  be  written  in  terms  of  integrals  over  phase  space 


V  V 
/  •  •/"\(0)\(t)e-^%^.  .  .d£^dr^.  .  .dr^ 

<V,(0)V^(t>=  ^^.^^ (5.7) 

'        '        -/SH 


e  ■  d£^...d£^dr^...dr^ 


-00    o 


Assumption  2  is  implicit,  since  the  integrals  are  performed  only  over 
center  of  mass  coordinates.  The  Hamiltonian  H  for  the  center  of  mass 
coordinates  is  given  by 


N   2 
H  =  \   TT^  +  U 


^2^^'  (5.8) 

1=1 
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and  if  Assumptions  3  and  4  are  satisfied,  we  may  write 


U  = 


N 

£  )        0 


L 


11 
a 


(5.9) 


Now,  if  we  define  reduced  quantities  as 


distance 


temperature 


time 


r 


T"  = 


a 

kT 
e 


,  1/2,  1/2 
=  te   /m  a 


(5.10) 
(5.11) 
(5.12) 


momentum 


-:^  ,,   111   111. 

3  =  p/(ni  '  £  '  ) 


(5.13) 


volume 


v"  =  v/cr 


(5.14) 


then 


-foo 


,v 


£ 
m 


^ 


p;(o)p:;(t")exp 


N   ':<1 
Pi 


<V  (0)V  (t)>  = 


-00    o 


V 


i=l 


ir 


■^     N^ 


exp 


-00    o 


N   -,v2 
2T* 


i=l 


N 


i<j 


d^i-.-d^dr^.-.d^ 


N 


—  N  0(r':.) 


d£^...d£^dr^...dr^ 


|g(T^v^t^"-) 

(5.15) 


where  g  is  a  universal  function  of  T  ,  V'  and  t  .   From  equation  (5.6) 


D  = 


^  ^g(T^v^t■"■)dt■" 

m     >^ 
o 


(5.16) 
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or  combining  equations  (5.15)  and  (5.16) 

D"  =  D"(t",v''')  (5.17) 


where         D"  =  ^  f^  (5.18) 


Similar  relationships  may  be  derived  for  shear  and  bulk 
viscosity  and  thermal  conductivity, 

r{'  =   Ti"(T",v'')  (5.19) 

/c*  =  k:"(T",v'')  (5.20) 

X''  =  xa*X')  (5.21) 

where  2 

n*  =^  (5.22) 

v  me 

/c"  =^  (5.23) 

/me 

X'^.^yY  (3.24) 

The  corresponding  states  principle  may  be  derived  under  less 
restrictive  conditions  than  those  used  above.   Although  it  is  then  of 
general  application  it  loses  the  attractive  simplicity  of  the  tijo- 
parameter  laws  expressed  by  equations  (5.17),  (5.19),  (5.20)  and 
(5.21).   Thus,  Assumption  1  is  readily  removed  by  using  quantum 
mechanical  expressions  in  place  of  integrals  over  phase  space.   This 
results  in  the  introduction  of  a  reduced  Planck  constant  in  the 
expressions  for  transport  properties. 
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where 


D"  =  D"(T",V",h')  (5.25) 

T)*  =    r]"(T",  V",  h")  (5.26) 

K*  =  k:"(T",  V",  h'")  (5.27) 

X*=   X"(t",  v",  h")  (5.28) 

u 

h"  = 


Donth  (77)  showed  that  the  corresponding  states  principle  could  be 
proved  for  thermodynamic  properties  without  assuming  pairwise  addi- 
tivity  of  potentials  (Assumption  3) .   Thus  in  place  of  Assumption  3 
he  assumed  that  the  total  potential  energy  of  the  system  might  be 
expressed  in  the  form 

U=e0(%  %....=!)  (5.29) 

where  0  is  a  universal  function  of  the  reduced  position  coordinates. 
Recently  Hakala  (75)  derived  a  more  general  form  of  the  corresponding 
states  principle  for  the  thermodynamic  properties  by  means  of  dimen- 
sional analysis.   The  effect  of  his  treatment  is  to  remove  Assumption 
4.   The  corresponding  states  principle  is  shown  to  apply  to  substances 
which  do  not  obey  a  pair  potential  of  the  form  given  by  equation  (5.1), 
provided  the  necessary  additional  reduced  parameters  are  included  in 
the  correlation. 

It  is  of  interest  to  note  that  the  above  corresponding  states 
relations  for  D  and  rj  may  be  derived  from  the  free  volume  equations. 
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The  assumptions  implied  are  in  this  case  less  clearly  defined.   The 
expression  for  self-diffusion  coefficient  derived  from  free  volimie 
theory  is  given  by  equation  (2.11)  as 


P ^     /^   ^v  \     /     ^o   \ 

^1/2  -  ^o^^P  V  ■  ^  J  ^""^  (  "  ^^^  )      ^^-^^^ 

Equation  (5.30)  can  be  expressed  in  terms  of  reduced  variables  as 


or 


where 


u       _-;<:/    V  \      /     o 


T  \   T         \   V  -  V 


o 


d''  =  D"(T",V")  (5.32) 


E 
^v  =  ^  =  eJ'(V*)  (see  Chapter  3) 

V 
Vq  =  "3  =  V^'(T")  (see  Chapter  3) 

B  m 
and  it  is  asstimed  that  B"  =  °  =  constant. 

Ok   ' 
in  this  chapter  we  consider  only  fluids  composed  of  monatomic 

molecules  (excluding  He)  so  that  Assumptions  1  and  2  may  be  regarded 
as  fully  satisfied.   The  inert  gas  molecules  have  been  found  to  obey 
the  principle  of  corresponding  states  for  thermodynamic  properties 
(20,74)  and  transport  properties  at  low  density  (78)  within  experi- 
mental accuracy.   Pairx^ise  additivity  of  potentials  seems  to  be  a 
good  assumption  up  to  reasonably  high  densities  for  the  molecules 
considered  in  this  chapter  (79-81).   Ross  and  Alder  (82)  have  pre- 
sented experimental  evidence  that  nonadditivity  becomes  a  significant 
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factor  for  these  fluids  only  at  very  high  densities,  of  about  4  times 
the  critical  value.   Assumption  4  concerning  the  form  of  the  pair 
potential  energy  is  only  approximately  satisfied  by  monatomic  mole- 
cules considered  here,  as  discussed  by  Smith  (83). 
5  .2 Previous  Work 

Previous  work  on  the  application  of  the  principle  to  trans- 
port coefficients  of  inert  gases  and  other  simple  polyatomic  molecules 
in  the  dense  gas  and  liquid  states  has  been  reported  by  Cini-Castagnoli, 
Pizzella  and  Ricci  (84),  Kerrisk,  Rogers  and  Hammel  (85),  Rogers  and 
Brickwedde  (86),  Boon  and  Thomaes  (87,88),  Boon,  Legros  and  Thomaes 
(64),  Tee,  Kuether,  Robinson  and  Stewart  (89),  Trappeniers,  Botzen, 
Van  Den  Berg  and  Van  Oosten  (90),  and  Trappeniers,  Botzen,  Ten  Seldam, 
Van  Den  Berg  and  Van  Oosten  (91). 

Cini-Castagnoli,  et_  al .  applied  the  principle  to  the  viscosity  and 
thermal  conductivity  of  saturated  liquid  Ar,  CO,  N„ ,  CH  and  0„ .  Correla- 
tions appeared  to  be  poor,  especially  for  the  reduced  viscosities. 

Kerrisk  et_   al.  and  Rogers  and  Brickwedde  examined  the  principle  of 

3    4 
corresponding  states  for  the  light  molecules  He  ,  He  ,  H» ,  D  ,  T  , 

etc.,  where  quantum  effects  are  important.   They  tested  the  principle 

for  the  saturated  liquid  viscosity  and  thermal  conductivity  of  these 

substances.   In  their  correlation  they  introduced  the  reduced  Planck 

constant  h  as  a  third  parameter  as  required  by  equations  (5.26)  and 

(5.28).   Boon  et  al .  tested  the  validity  of  the  corresponding  states 

principle  for  the  saturated  liquid  viscosity  of  Ar,  Kr,  Xe,  0_ ,  N  , 

C0_,  CH,  and  CD,.   They  found  that  the  inert  gases  obey  the  principle 


closely,  but  even  relatively  simple  polyatomic  molecules  such  as  0  , 

C'd/    and  CD,  were  found  to  deviate  markedly. 

Self-diffusion  coefficients  and  binary  diffusion  coefficients 

of  more  than  14  simple  molecules  in  the  dilute  gas,  dense  gas  and 

liquid  states  have  been  tested  by  Tee  _et  al.   They  defined  the  reduced 

diffusion  coefficient  d"   as 

AB 

T 
n   '   CAB 
c  D, 


"o  AB  V  P  ,^ 


AB    I ^  ^    ^  N  (5.33) 

"^cAB  V  m"  "*"  M~  j 
\  A    B  - 

where  c  here  represents  the  total  concentration,  and  obtained  from 

0 

the  equation  for  the  diffusion  coefficient  of  a  dilute  gas  the  relation 


D""  =  4.05  X  10"""  ^f^,„  (5.34) 

AB  ^(l,l)->' 

They  plotted  their  results  as  D"   against  P^  at  constant  T^  where 
^  ^  AB  R  R 

P  T 

P„  =  —  and  T„  =  —  .   However,  the  data  they  present  seem  to  be  too 

c  c 

scattered  to  justify  the  test. 

Trappeniers  e_t  al.  applied  the  corresponding  states  principle 
to  the  viscosity  of  inert  gases  up  to  high  densities  for  temperatures 
above  critical.   They  defined  the  reduced  viscosity  as 


Q*  =  2-  (5.35) 

o 


where 


°     16  TT  n^^'^^" 


and  2 

ri''=?i=.  (5.37) 
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Q"  was  then  related  to  p"  and  T" . 

5 ■  3 Test  of  Correspondence  Principle  for  Inert  Gases 

No  systematic  and  consistent  test  of  the  principle  for  all 
three  transport  properties  of  dense  fluids  over  a  wide  range  of  tempera- 
tures and  densities  seems  to  have  been  made.   In  particular  the  trans- 
port properties  of  fluids  at  very  high  pressures  (and  thus  high  densi- 
ties) have  not  been  examined.   The  behavior  at  high  densities  is  of 
particular  interest  since  it  is  under  these  conditions  that  Assumption 
3  may  become  invalid.   Also,  errors  arising  from  differences  in  the 
repulsive  portion  of  the  intermolecular  potential  for  the  molecules 
considered  (Assumption  4)  may  become  apparent  at  high  density. 

The  essential  simplifying  factor  in  the  case  of  monatomic  (as 
opposed  to  polyatomic)  molecules  is  that  Assumption  2  is  fully 
satisfied.   A  comparison  of  the  correspondence  behavior  for  monatomic 
and  polyatomic  molecules  may  therefore  throw  light  on  the  validity  of 
this  assumption.   Calculations  are  presented  in  this  chapter  to  test 
the  correspondence  principle  for  all  three  transport  coefficients  over 
the  full  range  of  temperatures  and  pressures  for  which  data  are 
available.   The  extension  to  polyatomic  molecules  is  discussed  in  the 
following  chapter. 

Reduction  of  the  transport  coefficients  are  performed  using 
Lennard- Jones  [6,12]  potential  parameters  rather  than  critical  constants 
because  the  former  have  more  theoretical  significance.   The  use  of 
critical  constants  is  justified  when  the  group  of  substances  accurately 
obey  the  same  potential  function,  Assumption  4  (50),   In  such  a  case 
it  would  in  principle  be  immaterial  whether  potential  function  or 
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critical  parameters  were  used.   However,  in  practice  the  group  of 
molecules  do  not  all  accurately  obey  Assumption  4,  and  the  potential 
parameters  are  force-fitted  to  some  semi-empirical  equation,  in  this 
case  the  Lennard- Jones  [6,12]  potential  function.   Under  such  condi- 
tions the  use  of  potential  parameters  rather  than  critical  constants 
may  improve  the  correlation  obtained.   In  this  worl<,  reduction  with 
both  critical  constants  and  potential  parameters  (smoothed  in  the  way 
described  below)  was  tried.   The  use  of  potential  parameters  was  found 
to  give  consistently  better  results. 

The  potential  parameters  reported  in  the  literature  by  differ- 
ent workers  show  substantial  variations.   Thus  for  krypton  values  of 
£ ,  the  characteristic  energy,  reported  by  different  authors  vary  by 
as  much  as  14%.   It  is  therefore  desirable  to  find  some  means  of 
averaging  these  parameters.   Moreover,  Reed  and  McKinley  (92)  have 
shown  that  polyatomic  molecules  do  not  all  obey  the  simple  two-para- 
meter Lennard- Jones  [6,12]  potential  energy  function.   Since  Assumption 
4  requires  that  all  molecules  should  obey  the  same  potential  energy 
function  in  order  to  arrive  at  the  same  corresponding  states,  an 
attempt  was  made  to  force  the  polyatomic  molecules  to  fit  a  two- 
parameter  potential  energy  function  by  adopting  the  method  used  by 
Tee,  Gotoh  and  Stewart  (93).   According  to  the  Pitzer  compressibility 
correlation  (94)  the  volumetric  behavior  of  a  fluid,  and  thus  the 

intermolecular  forces  present,  are  characterized  by  P  ,  T  and  the 

c   c 

acentric  factorco .   Thus  the  intermolecular  force  constants  e  and  a 
should  depend  on  these  three  quantities,  so  that  one  obtains  the 
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dimensionless  relations 


^  =  fCcu)  (5.38) 

c 


P   N  1/3 

^    )         a  =  f  (CO)  (5.39) 

c  ^ 


where  f  and  f  are  universal  functions.   In  Figures  5.1  and  5.2 

values  of  (  —2—   ^      and  :; for  twenty  different  molecules  are 

V  kT    y  °       kT 

^    c  ■/  c     , 

plotted  against  their  corresponding  Pitzer  factor  co.   The  best  straight 


lines  X'/ere  drawn  through  them  by  the  method  of  least  squares.   Poten- 
tial parameters  used  in  these  figures  were  values  calculated  from  gas 
viscosity  data,  and  were  taken  from  Hirschfelder ,  Curtiss  and  Bird 
(20)  and  Tee,  Gotoh  and  Stewart  (93).   Where  several  values  fore  and  a 

were  reported  for  a  given  molecule  the  mean  values  are  plotted  in 
Figures  5.1  and  5.2,  and  the  maximum  deviation  is  indicated.   The 
calculations  presented  in  this  chapter  and  the  next  make  use  of  smoothed 
potential  parameters  obtained  from  these  two  figures.   Values  of  the 
potential  parameters  used  are  shown  in  Tables  5.1  and  6.1.   Values  for 
molecules  other  than  those  tabulated  may  be  calculated  from  the 
equations  for  the  best  lines  shown  in  Figures  5.1  and  5.2  which  are 

-^  =  0.7932  -  0.672300  (5.40) 

kT 
c 

^  Y/2^=  0.4583  +  0.1213a)  (5.41) 

kT   ] 
c  ^ 

Figures  5.3,  5.4  and  5.5  are  reduced  plots  of  saturated  liquid 
viscosity,  thermal  conductivity  and  self-diffusion  coefficient  as 
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functions  of  temperature  for  the  inert  gases  (argon,  krypton,  xenon 
and  neon) .   For  pure  saturated  liquids  there  is  only  one  degree  of 
freedom,  and  the  reduced  transport  coefficients  depend  only  on  reduced 
temperature.   In  the  saturated  liquid  region,  the  inert  gases  appear 
to  obey  the  corresponding  states  principle  for  the  above  transport 
properties  reasonably  well  with  the  exception  of  the  thermal  conducti- 
vity data  of  neon.   The  discrepancy  observed  in  Figure  5.4  for  neon 
may  be  due  to  quantum  effects,  since  neon  is  the  second  lightest 
molecule  among  the  inert  gases.   Experimental  data  for  bulk  viscosity 
are  scarce.   The  only  bulk  viscosity  measurements  made  seem  to  be 
those  reported  by  a  few  workers  for  argon  (14,95,96).   Figure  5.6 
is  a  reduced  plot  of  bulk  viscosity  of  saturated  liquid  argon 
(extrapolated  data)  as  a  function  of  temperature.   Great  experimental 
difficulties  are  involved  in  bulk  viscosity  measurements,  and  the  data 
are  very  scattered.   Even  among  the  data  of  the  same  worker  the 
average  percent  deviation  is  as  much  as  157<,. 

In  the  derivation  of  the  correspondence  principle  for  trans- 
port properties  shown  above,  reduced  viscosity  is  found  to  be  a 
function  of  reduced  temperature  and  reduced  molal  vol;ame.   Thus 
according  to  equation  (5.19)  T]'   =  r)^'(T'%V").  -However  it  is  more 
convenient  in  practice  to  plot  the  transport  properties  as  a  function 
of  temperature  and  pressure.   Since  reduced  molal  volume  may  be  shown 
to  be  related  to  reduced  pressure  through  the  equation  of  state,  one 
may  express  t}"  as  r)"(T",P"). 

Reduced  viscosities  of  inert  gases  plotted  as  a  function  of 
reduced  temperature  at  constant  reduced  pressures  are  shown  in 
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Figure  5.7.   The  corresponding  states  principle  appears  to  apply 
closely  to  viscosity  of  inert  gases  over  the  entire  range  of  tempera- 
ture and  pressure  for  which  experimental  data  are  available.   Data 
plotted  in  Figure  5.7  cover  densities  from  the  dilute  gas  phase  to 
the  dense  gas  and  liquid  phase,  and  correspond  to  pressures  from  zero 
to  over  2000  atmospheres.   A  similar  plot  for  the  thermal  conductivity 
of  inert  gases  is  shown  in  Figure  5.8.   In  this  plot  the  correspondence 
principle  is  found  to  apply  well  for  argon,  krypton  and  xenon  over 
the  entire  temperature  and  pressure  range  for  which  data  are  avail- 
able.  The  neon  data  in  Figure  5.8  obey  the  simple  correspondence 
principle  at  high  temperatures,  but  not  at  low  temperatures.   This 
behavior  for  neon  suggests  that  quantum  effects  are  important  for 
this  fluid  at  low  temperatures. 

The  reduced  viscosity  and  reduced  thermal  conductivity  data 
are  found  to  fit  equations  of  the  form 

In-rj''  =  a^  +  a^T""^  +  a^T*"^  +  a^f''"^  +  a^T*"^     (5.42) 


,* 


■"'-1   ,   K   t'''-2   ,   -L   rp"V-3       ,    rp*-^ 


and  InX  =  ^^  +  b2T""-^  +  b^T  "^  +  b^T  -^  +  b^T""^     (5.43) 

at  constant  reduced  pressure  or  along  the  saturation  curve.   The 
reduced  self-diffusion  coefficient,  on  the  other  hand,  is  found  to 
fit  an  equation  of  the  form 

InD*  =  d  +  d  t""^  (5.44) 

The  principal  parameters  used  to  test  the  correspondence 
principle  for  monatomic  molecules  are  shown  in  Table  5.1.   The 
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potential  parameters  included  in  this  table  are  obtained  from 
Figures  5.1  and  5,2,   The  coefficients  of  equation  (5,42)  evaluated 
by  computer  for  saturated  liquid  viscosity  of  monatomic  molecules 
(Ar,  Kr,  Xe,  Ne)  together  with  the  average  error  between  the  equation 
and  experiment  are  shown  in  Table  5,2,   The  average  percentage  errors 
betivfeen  prediction  and  experiment  for  the  inert  gases  are  within  kX, 
and  this  is  of  about  the  same  order  as  the  experimental  accuracy. 
Tables  5,3  and  5.4  exhibit  the  coefficients  of  equations  (5.43)  and 
(5.44)  respectively  for  monatomic  fluids.   The  maximum  disagreement 
between  reported  thermal  conductivity  data  of  different  workers  for 
monatomic  molecules  is  estimated  to  be  approximately  1°L.      Equation 
(5,43)  X'/hen  used  with  the  values  of  coefficients  shov/n  in  Table  5,3 
gives  an  average  percentage  error  well  within  47o  (thermal  conductivity 
data  of  neon  were  not  used  to  fit  the  coefficients) ,   Disagreements 
between  self-diffusion  coefficient  data  of  different  workers  are 
especially  large.   Thus  the  seif-dif fusivity  of  xenon  reported  by 
two  v7orkers  differ  by  as  much  as  100%  (18,57).   Equation  (5.44) 
for  self-dif fusivity  of  monatomic  molecules  using  the  values  of 
coefficients  in  Table  5.4  yields  an  average  percentage  error  of  about 
5%.   Only  the  data  of  Naghizadeh  and  Rice  (18)  were  used  in  obtaining 
these  coefficients. 

Table  5.5  lists  the  coefficients  of  the  viscosity  equation 
together  v;ith  percentage  errors  for  the  viscosity  isobars  shown  in 
Figure  5.7.   The  average  errors  are  largest  at  the  highest  reduced 
pressures,  v/hich  correspond  to  pressures  up  to  approximately  2000 
atmospheres.   These  larger  errors  may  be  due  to  a  failure  of 
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Assijmption  3  at  these  high  pressures,  or  to  experimental  errors.   The 
coefficients  of  the  ther.:ial  conductivity  equations  for  the  isobars 
shCT^7n  in  Figure  5.8  and  the  average  percentage  error  between  experi- 
ment and  the  equations  are  given  in  Table  5.6.   The  errors  for  all  the 
equations  are  within  experimental  error.   Sources  of  viscosity, 
thermal  conductivity  and  self-diffusion  coefficient  data  are  also 
included  in  Tables  5.2  -  5.6. 

5  .4 Summary 

The  good  agreement  of  the  experimental  data  with  the  principle 
of  corresponding  states  leads  one  to  conclude  that  the  4  assumptions 
on  which  the  molecular  derivation  is  based  are  good  approximations 
for  monatomic  fluids.   An  exception  is  neon,  for  which  quantum  effects 
must  be  taken  into  account  at  low  temperatures.   It  is  of  particular 
interest  that  Assumption  3  concerning  pairvise  additivity  of  the 
intermolecular  potential  does  not  lead  to  significant  errors  even  at 
the  highest  densities  for  which  measurements  are  available.   Similarly, 
Assumption  4  does  not  give  rise  to  significant  errors  for  these 
fluids  in  spite  of  known  differences  in  the  potential  energy  form 
(S3). 
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CHAPTER  6 

PRINCIPLE  OF  CORRESPONDING  STATES  FOR  POLYATOMIC 
NONPOLAR  FLUIDS 


It  has  long  been  realized  that  the  simple  correspondence 
principle  does  not  apply  accurately  to  polyatomic  fluids  in  the  con- 
densed phase.   In  the  previous  chapter  it  has  been  shown  that  monato- 
mic  molecules  obey  the  principle  of  corresponding  states  over  the 
entire  range  of  temperature  and  pressure  conditions.   In  this  chapter 
a  careful  re-examination  of  the  four  basic  assumptions  used  in  deri- 
ving the  principle  is  first  made.   A  simple  model  of  the  molecular 
behavior  of  a  polyatomic  liquid  is  proposed,  and  provides  the  basis 
for  a  rederivation  of  the  correspondence  principle. 
6  . 1 Problems  in  Polyatomic  Fluids 

Polyatomic  fluids  should  obey  Assumption  1  closely  since  their 
molecular  weights  are  usually  large  and  thus  quantiom  effects  are 
negligible.   Assumption  3  concerning  pairwise  additivity  of  potentials 
does  not  seem  to  be  a  major  source  of  error  for  either  monatomic  or 
polyatomic  fluids  (81) .   Therefore  it  appears  likely  that  the  deviation 
of  polyatomic  fluids  from  the  correspondence  principle  is  due  to  the 
fact  that  they  do  not  satisfy  either  or  both   .  Assumptions  2  and  4. 

T\^70-parameter  potential  energy  functions,  such  as  the  Lennard- 
Jones  [6,12]  equation,  do  not  accurately  account  for  properties  of 
polyatomic  fluids  (92).   In  fact  three-parameter  energy  functions,  such 
as  the  Lennard- Jones  [6,n]  or  the  Kihara  potential  give  better  results 
for  these  molecules.   Thus  polyatom.ic  fluids  in  general  do  not 
satisfy  Assumption  4  closely.   Failure  of  polyatomic  fluids  to  satisfy 
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Assimiption  4  may  be  corrected  by  introducing  a  suitable  third  para- 
meter.  Alternatively  a  tt-.'o- parameter  potential  may  be  used,  taking 
suitable  potential  parameters  in  the  way  suggested  by  Tee,  Gotoh  and 
Stewart  (93). 

Recently  Preston,  Chapman  and  Prausnitz  (112)  examined  the 
corresponding  states  principle  for  the  transport  coefficients  of  nine 
monatomic  and  polyatomic  fluids  in  the  saturated  liquid  state.   They 
proposed  that  potential  parameters  of  substances  in  the  liquid  state 
may  take  different  values  from  those  in  the  gaseous  state.   Liquid 
phase  potential  parameters  x^7ere  obtained  from  the  saturated  liquid 
viscosity  and  thermal  conductivity  data  by  assimiing  that  the  simple 
correspondence  principl.^  applied.   They  then  used  these  parameters  to 
test  the  correspondence  principle  for  the  saturated  liquid  viscosity, 
thermal  conductivity  and  self -diffusion  coefficient  of  the  nine 
molecules.   This  procedure  gave  reasonably  good  correlations  for  liquid 
viscosity  and  thermal  conductivity,  with  root-mean- square  deviations 
of  7.47o  and  5.4%  respectively.   The  correlation  for  self-diffusion 
coefficient  was  comparatively  poorer,  showing  a  root-mean- square  devia- 
tion of  157c.   Their  treatment  seems  open  to  criticism  on  several 
counts.   Firstly,  by  force-fitting  the  parameters  using  liquid  phase 
data  the  validity  of  their  test  is  doubtful.   In  particular,  certain 
facets  of  the  behavior  of  polyatomic  molecules  may  be  hidden  in  this 
way;  thus  hindered  rotation  in  the  liquid  state  may  be  important,  as 
discussed  below.  Further,  because  separate  sets  of  potential  parameters 
are  iiecessary  for  the  gas  and  liquid  phases,  problems  will  arise  in 
choosing  the  most  suitable  parameters  for  the  intermediate  densities, 
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and  it  will  not  be  possible  to  predict  transport  properties  when  no 
experimental  liquid  phase  data  are  available. 

Little  attention  seenis  to  have  been  given  previously  to  the 
possible  failure  for  liquids  of  Assumption  2  concerning  the  aensity 
independence  of  intramolecular  degrees  of  freedom.   Because  vibration- 
al quanta  and  bond  energies  are  large  one  would  not  expect  inter- 
molecular  forces  to  have  much  effect  on  the  vibrational  energy  of  a 
polyatomic  molecule.   However,  rotational  energies  are  smaller  and  it 
seems  quite  likely  that  free  rotation  of  polyatomic  molecules  may  not 
be  possible  at  high  densities. 

Figure  6.1  shows  a  plot  of  r)"  at  P"  =  0  for  a  variety  of 
monatomic  and  polyatomic  molecules  in  both  the  gas  and  liquid  states. 
Reduction  of  the  dilute  gas  and  liquid  viscosity  coefficients  is 
performed  using  smoothed  molecular  parameters  determined  from  gas 
viscosity  data  (see  Figures  5.1  and  5.2).   In  the  dilute  gas  phase 
both  monatomic  and  polyatomic  molecules  obey  the  correspondence 
principle  very  well.   In  the  liquid  phase,  however,  only  the  monatomic 
molecules  obey  the  principle.   Deviations  for  the  polyatomic  molecules 
increase  as  the  density  increases.   Were  these  deviations  due  only  to 
incorrect  potential  parameters,  as  suggested  by  Preston  et  al.  (112), 
one  would  expect  that  discrepancies  might  occur  for  monatomic  molecules 
in  the  liquid  state  also.   The  fact  that  gas-phase  parameters  work 
well  for  the  latter  molecules  in  the  liquid  states,  together  with  the 
excellent  agreement  found  for  polyatomic  fluids  at  low  density, 
strongly  suggests  a  failure  of  Assumption  2  in  the  liquid  phase.   It 
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seems  likely,  therefore,  that  hindered  rotation  has  a  significant 
effect  on  the  transport  coefficients  of  polyatomic  fluids  at  high 
densities,  and  that  this  effect  causes  increasingly  large  deviations  as 
density  rises.   Davis  and  Matheson  (113)  have  presented  evidence  to 
show  that  rotational  degrees  of  freedom  are  restricted  at  high  densi- 
ties. 

Figure  6.2  shows  reduced  viscosity  for  20  saturated  liquids. 
The  same  trends  as  in  Figure  6.1  are  seen;  deviations  from  the  curve 
for  monatomic  liquids  are  greatest  for  molecules  having  larger  Pitzer 
factors . 

6.2 Derivation  of  Corresponding  States  Principle  for  Polyatomic 

Molecules 

In  the  gas  phase  molecular  separations  are  usually  large  and 
intermolecular  forces  are  small.   A  molecule  can  therefore  rotate 
freely  about  its  body  axes.   However,  as  the  density  increases  inter- 
molecular forces  become  more  important,  and  if  these  forces  are  non- 
central  the  molecules  exert  a  torque  on  one  another.   Polyatomic 
molecules  thus  rotate  in  a  potential  energy  field  which  is  a  function 
of  the  relative  orientations  and  center  of  mass  positions  of  the 
molecules;  in  the  case  of  very  high  densities  the  rotational  motion 
might  be  reduced  to  that  of  angular  vibrations  about  the  body  axes. 

For  dense  fluids  the  rotational  degrees  of  freedom  thus  become 
dependent  on  the  center  of  mass  coordinates,  and  cannot  be  factored 
out  from  the  phase  integral  as  in  the  case  of  dilute  gases.   The  time 
correlation  function  instead  should  be  written  as 
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•H/kT , 

X      X 


V^.(0)V^.(t)e     ^d£^  •■•dV^e  •••d0^ 


1 


<V  (0)V  (t)>  = 
^  X     X    ^ 


../  e      d£^  •••dr^d£g  .-.d^^  (^-l) 


1 


where   H(p   ...r  ..pn  — Q)    =  K„  (p   ...)  +K   (p^  ...) 
r^       -N    ^2.   ~^      CM^-'=^r,    '^     rot^-^&    '^ 

+  U(r^...0^...) 

The  phase  integral  in  this  form  will  be  difficult  to  handle, 
since  we  have  little  knowledge  as  to  the  form  of  the  potential  energy 
term.   Although  the  kinetic  energies  of  translation  (Kpv)  ^^'^   rotation 
(K   )  are  separable,  the  potential  energy  term  is  not.   Evaluation 
of  the  integrals  therefore  requires  that  a  suitable  expression  be 
obtained  for  U(_r,_0)  .   The  part  of  the  Hamiltonian  arising  from  intra- 
molecular vibrations  is  here  assumed  to  be  separable  and  independent 
of  density  for  the  reasons  discussed  above. 

To  account  in  an  approximate  way  for  the  coupling  between 
translational  and  rotational  degrees  of  freedom,  a  highly  simplified 
one-dimensional  model  is  proposed.   Consider  a  harmonic  oscillator 
cell  model  similar  to  the  type  used  by  Prigogine  (50)  .   In  this  model 
each  molecule  is  assijmed  to  spend  a  large  fraction  of  its  time  oscil- 
lating in  a  cell  or  cage.   Since  the  oscillator  is  assumed  harmonic, 
the  restoring  force  acting  on  a  ir.olecule  at  any  time  is  proportional 
to  the  displacement  of  the  molecule  from  its  equilibrium  position  along 
each  of  its  orthogonal  coordinates.   The  nearest  neighbor  molecules 
are  assumed  in  an  average  situation  to  be  on  their  lattice  sites,  and 
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are  treated  as  smeared  over  a  sphere  passing  through  their  lattice 
positions.   For  the  present  a  one-dimensional  motion  is  assumed. 
However,  it  can  be  readily  extended  "o  the  three-dimensional  case  by 
assuming  an  isotropic  system  in  which  the  motion  along  any  of  the  3 
axes  is  the  same. 


Figure  6.3.  A  One-Dimensional  Oscillator  Model. 


Consider  a  molecule  oscillating  linearly  along  the  x-axis  in 
a  cell  (Figure  6.3).  For  convenience  the  molecule  is  represented  as  a 
rigid  diatomic  one,  but  similar  considerations  will  apply  in  general. 
Because  of  the  linear  oscillation,  the  molecule  rotates  in  an  oscil- 
lating potential  field,  and  this  gives  rise  to  hindered  rotation. 
The  rotational  motion  of  the  molecule  about  the  z-axis  (perpendicular 
to  the  plane  of  the  paper)  may  be  considered  to  consist  of  tvjo  parts: 
a  constant  velocity  rotation  plus  a  harmonic  angular  oscillation 
superimposed  on  it.   The  harmonic  angular  oscillation  is  assumed  to  be 
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in  phase  with  the  linear  oscillator.   The  constant  velocity  rotation 

is  the  value  that  would  be  observed  if  no  linear  oscillation  occurred, 

that  is  when  the  molecule  remained  at  0.  The  angular  coordinates  can 
be  written  as 


e   (x,t)  =  (D  (t)  +  ¥  (x)  (6.2) 

XXX 

^e^'^^  =  P(i)   "^  P^  (^^  (6.3) 

X  XX 


where  0  (t)  and  p-  refer  to  the  constant  velocity  rotational  motion 

X 

at  0,  while  ^  (x)  and  p   (x)  account  for  the  superimposed  oscillation 

X 

or  angular  perturbation. 


Pg  =P0  +  P^ 

X     X   ^  X 


rotating  reference 
axis 


Figure  6.4.   Superimposed  Rotational  Motions 
of  a  Molecule. 


The  coordinates  are  shown  in  Figure  6.4  for  an  instant  in  time  when 

the  molecular  axis  lags  behind  the  rotating  reference  axis.   The 

reference  axis  rotates  with  the  constant  velocity  C>  ,  and  viewed 

x 
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relative  to  this  axis  the  molecule  performs  angular  vibrations. 

Because  the  t^^ro  oscillations  are  in  phase  the  quantities  ^  and  p„, 

X      Y 

X 

will  depend  only  on  the  linear  displacement  x. 

Using  equations  (6.2)  and  (6.3)  it  is  possible  to  transform 

from  the  phase  space  coordinates  p  ,  0  used  in  equation  (6.1)  to 

a     — 

the  new  angular  coordinates  p ,,  0 ,  so  that 

-=-q)  — 


<v^(o)v..(t)>  = 


|v^(0)V^(t)e-^^/^^  d^     dr^dp     dO^ 


^.Te-^/^^dp^    dr^d£^    dO^ 
^  11 


(6.4) 


where   H  =  H(p^,  1,    V^^ 

and  it  is  noted  that  for  coordinate  transformation  in  phase  space  the 
Jacobian  is  unity  (114) .   In  practice  the  above  transformation  requires 
knowledge  of  the  functions  Y   (x)  and  p^  (x) .   Expressions  for  these 

X 

quantities  are  obtained  below. 

For  small  oscillations  the  total  potential  energy  can  be 
expressed  in  terms  of  a  Taylor's  series  expansion  as 


\    y    x=0    \   X  /  x=0         \  dx   x=0 

Y  =0  ^  =0  Y  =0 

XX  X 

\  ^^^^x  y  x=o  ^    ^\b^^  J 

^  =0  ^   x=0 

^  w  =0 

X 
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Since  U(x,Y  )  is  a  minimum  at  x  =  0  and  ^  =  0 

X  X 


f^^ 


8u 


V  ox  /   ^   V  o"?  . 

„x=0    V   X  /  x= 

X  ^  X 


=  0 


Denoting 


(6.6) 


Neglecting  terms  higher  than  second  order,  \-je   have 

^  ox  /  x=0      \     X  /  x=0 
^  2  (  ~  )    \  (6.7) 

Vox  ^^^^Q 

^2=  (Ur}  (6-9) 

x=0 

^22  =  (  7~2  )  (6.10) 

x=0 
equation  (6.7)  can  be  written  as 

U(x,^if^)  =  U(0,0)  +  I  k^^x^  +  k^2^Y^  +  I  k22T^  (6.11) 

For  small  oscillations  the  total  kinetic  energy  may  be  written  as 

1-2    1   -^ 
K  =  -  mx^  +  -^  lY  (6.12) 

2       2   X  ' 

■    _    dx  *  dY 

where  x   -  -^  and   ^J/^  =  —  ,    and    the      Lagrangian  of   the   system  becomes 


L  =  K(x,Y   )    -   U(x,Y   ) 


1      .2 

—  mx 


+  ^  I^^  -    [U(0,0)   +1  k^^x^   +  k^2^Y^  (6.13) 


+  I  ^22^] 
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The  Xa'-i^rangian  equation  of  motion  with  respect  to  x  is  given 
by 

dt  \  ox  /    ox 
So  that  from  equation  (6.13) 

mii  +  k^  'X  +  k.-^   =  0  (6.15) 

Similarly  the  Lagrangian  equation  of  motion  with  respect  to  ^  is 

dt  U\^    V  a\^  (6.16) 


or 


Let 


I^  +  k„.^  +  k..x  =  0  (6.17) 

x     22  X     12 


X  =  A'cosOO  t  +  A"sinXi,t  (6.18) 

^"^^  ^     =  B'cosco.t  +  B"siilO,t  ^^-^^^ 

X  1  1 

be  the  solutions  to  equations  (5.15)  and  (6.17).   A  complete  solution 
of  equations  (6.15)  and  (6.17)  is  shown  in  Appendix  3.   Differentia- 
ting equations  (6.18)  and  (6.19)  twice  with  respect  to  time  yields 

X*  =  -CO,  [A'cosoo^t  +  A"sina),t] 
1       i         1 

=  -CD^x  (6.20) 

and  2 

¥  =  -CUUB'cosoo^t  +  B"sim),t] 
X     1       1  1 

=  -a)?7  (6.21) 

1  X 

Substituting  equations  (6.20)  and  (6.21)  into  equations  (6.15)  and 
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(6.17) 


?mx  +  k,.x  +  !<,,•?     =0  (6.22) 


-cc^r?     +  k-.Y     +k.-x  =  0  (6.23) 

i     X  22  X  12 


2 


Eliniinating  CD     and   rearranging  we  have 


9 

IT-  ^  2 


mx 


c 

The  potential  energy  U(x,T  )  now  becomes 


U(x,-^)  =  U(0,0)  +1  [k^^x^  4-  k^2-\]  +i  ^^22^^  +  k^2^^?J 

mx  (6.25) 


Let  us  consider  at  any  time  t  =  t   x  =  0  and  Y   -  0,  then 

0  =  A'cosa>[_t^  +  A"sina>[_t^  (6.26) 

^  =   -cotCO^t  (6.27) 

A         1  o 


and 


0  =  B'cosoo.  t  +  B"sinait  (6.28) 


or 


therefore 


1  o 


p-  =  -cotco^t^ 


(6.29) 


A^  ^  b;^  (6.30) 

A'    B' 


Combining  equations  (6.18)  and  (6.19)  we  have 

B'cosCD-jt   +  B"sinJa  t 
•a      = X 

X       A'cosco,t  +  A"sina)  t 
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B" 
B'[cosco  t   +-7  sincu  t] 
i B 1 

A'Lcosco^t  +  YT  Sim),  t] 
1  A  1 


^  (6.31) 


where  „  _  B' 

""  "  A^  (6.32) 


Equation  (6.25)  can  be  rewritten  as 

U(x,^  )  =  U(0,0)  +  [  ^=^  +  1]  -i  [k^.x^  +  k  .x^  ]  (6.33) 

>-  m       2    11      12   X 

Differentiating  equation  (6.31)  with  respect  to  time  we  get 


?^  =  Sx  (6.34) 

The  total  kinetic  energy  is  given  by 

K  =^  mi^  +-^  I^^  (6.35) 

2       2   ^x 

and  making  use  of  equation  (6.32)  we  have 

^    1   .2  _^  1  ^r^.T2 

K  =  -T  mx  +  -r  ILcxJ 

„2 

=  1  mx^[  ^^-^   1]  (6.36) 

Now  we  define 

c  =  [i^  +  1]  (6.37) 

m 
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The  moment  of  inertia  appearing  in  the  above  equations  is  that 
for  the  molecule  rotating  about  an  axis  fixed  in  space,  rather  than 
about  a  body  axis.   It  may  be  regarded  as  a  time-average  of  the 
principal  m.oments  of  inertia  for  the  molecule.   In  the  limit  where 
hindered  rotation  is  absent,  as  for  dilute  gases  and  monatomic  molecules 
at  all  densities,  B'  and  B"  are  zero  and  so  is  S,  so  that  c  reduces 
to  unity. 

The  Hamiltonian  for  the  system  of  N  molecules  in  one  dimension 
can  be  written  as 

^  cp2  ^  4  i 

=  ^_lx^eU(x^...^)+  ^  _^  ^U(O) 

i=l  ■   i=i  (6.38) 


where 


N  N 


U(x^...x^)  =\    V  i\A^\2-^\^-\    -\    fkii  +-^2^  X  4 
i=l  i=l 


and  U(0)  is  the  potential  energy  when  all  molecules  are  at  their  cell 

\~'   2 
centers.   In  the  above  equation   >   p.  /2m  represents  the  kinetic 

energy  due  to  oscillations  in  the  x-direction,  and  U(x  . . -x  )  is  the 

potential  energy  associated  with  the  linear  oscillations.   The  term  c 

thus  accounts  for  energies  associated  with  the  angular  vibrations.   The 

last  two  terms  in  equation  (6.38)  corresponding  to  the  cj)  coordinate 

are  now  separable,  and  the  center  of  mass  Hamiltonian  for  the  one 

dimensional  case  is 

N    2 

\   =    l_         1^     +CU(X^...X^)  ^g_3g) 

i=l 
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Comparison  of  this  equation  with  equation  (5.8)  for  the  center  of  mass 
Hamiltonian  when  there  is  free  rotation  shows  that  the  two  expres- 
sions are  of  the  same  form.   In  fact,  as  far  as  the  equations  of 
classical  mechanics  are  concerned,  we  may  treat  the  hindered  rotation 
case  in  the  same  way  as  free  rotation,  by  replacing  the  usual  kinetic 

and  potential  energy  terms  by  those  in  equation  (6.39),  that  is  by 

N 

)   cp   /2m  and  cU  respectively.   The  new  kinetic  energy  term  implies 
A    ix 

1^  1/2 

that  momentum  p  should  be  replaced  by  c   p. 

Assuming  the  oscillation  is  isotropic,  so  that  the  motions 

of  the  molecule  along  all  the  3  orthogonal  coordinates  are  equivalent, 

the  time  correlation  function  can  now  be  written  as 


f.  .  .  rcV^(0)V^(t)e"^CM/'-T  ^^^  _  _  -dp^dr^.  .  .dr^ 


<y^^^K^^^>  =  S-f 


J...  I   e  d£^...d£^.dr^...dr^ 

-00     o 

3N   2 

Zcp. 
^  +  cU(r^...V 
i=l 
Now  define  the  reduced  temperature  T  as 


(6.40) 


T"  =  —  (6.41) 

ce 


where  c  is  given  by  equation  (6.37).   Writing  equation  (6.40)  in  terms 
of  reduced  variables,  yields 


^The  derivation  for  the  three-dimensional  case  is  considered 
in  more  detail  in  Appendix  4. 
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<V_(0)V  (t^ 


-TOO         V 

'I    I-f 

-00          o 

p;:(o)p;(t*)exp 

N       .v2 

-y '' 
4  ^^^ 

+00        V" 

/ . . . /   exp 
J       J 

-00          o 

N          ,V9 

L         L    2T* 

1=1 

N 

-^y  <3>i^£t---^% 


f"  A-       'J    J 

1^1 


•N 


=—  g(T",V",t) 

1    O    -"-   1   -"-     ■'<  ""  (6-^2) 

^  ..  Kj  .    -^ 

where  V",  p",  r"  and  t"  are  as  defined  in  Chapter  5.   The  self- 
diffusion  coefficient  is  given  by 


1/2   ~    , 
D  =  ^^-YJT      /  <\(0)\(t'"')>  dt"  (6.43) 


or 


^-  =   D"(T",V")  (6.44) 


where  D"  is  given  by 

D  =-—./-  (6.45) 

cO\/  £ 

Similarly  it  can  be  shown  that 

r^'  -Tf'C^'X)  (6.46) 


k:  -  /c  (T  ,V  ; 


(6.47) 


A.  =  X  (T  ,V  )  (6.48) 


where  now  the  reduced  shear  viscosity,  bulk  viscosity  and  thermal 
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conductivity   are   given  by 


'n     =  ~^ZZ  (6.49) 


/me 


2 
k:     =  (6.50) 


c 


/me 


'^ 


For  dilute  gases  and  monatomic  molecules  c  reduces  to  unity 
and  therefore  equations  (6.45),  (6.49),  (6.50)  and  (6.51)  all  reduce 
to  the  forms  obtained  by  the  simple  correspondence  principle  derivation 
in  equations  (5.18),  (5.22),  (5.23)  and  (5.24).   Thus  the  above  is  a 
general  derivation  of  the  correspondence  principle  which  should  apply 
to  dilute  and  dense  phases,  and  to  both  monatomic  and  polyatomic 
nonpolar  fluids. 
6 . 3 Test  of  Correspondence  Principle  for  Polyatomic  Molecules 

Several  workers  (87,88)  have  pointed  out  that  the  correspond- 
ing states  principle  in  its  simple  form  does  not  apply  to  transport 
properties  of  even  relatively  simple  polyatomic  molecules,  such  as 
nitrogen  and  carbon  monoxide,  in  the  dense  phase.   Some  investigators 
have  adopted  the  pessimistic  viewpoint  that  the  principle  is  inapplicable 
for  such  fluids.   In  this  section  the  modified  corresponding  states 
principle  derived  above  is  to  be  tested  for  the  viscosity,  thermal 
conductivity  and  self-dif fusivity  of  polyatomic  molecules  over  a  wide 
range  of  temperature  and  density  conditions. 

To  make  a  test  of  the  proposed  equations,  reasonably  accurate 
values  of  the  characteristic  factor  c  of  different  fluids  must  be 
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determined.   This  factor  can  be  evaluated  from  one  experimental  value 
for  any  one  of  the  transport  coefficients.   Saturated  liquid  viscosi- 
ties were  used  to  determine  c,  as  they  seem  to  be  the  most  accurately 
measured  of  the  transport  coefficients,  and  are  most  readily  available 
in  the  literature.   The  factor  c  of  25  monatomic  and  polyatomic  mole- 
cules thus  evaluated  is  shown  in  Table  6.1.   A  complete  list  of  the 
smoothed  potential  parameters  evaluated  from  Figures  5.1  and  5.2  is 
also  included  in  Table  6.1.   Table  6.2  gives  the  sources  of  transport 
property  data  used  here  to  test  the  correspondence  principle. 

The  viscosity  data  for  the  low  boiling  molecules  are  estimated 
to  have  a  maximum  disagreement  bet^.een  workers  of  approximately  81; 
an  exception  is  oxygen  for  which  the  maximum  disagreement  is  20%. 
Saturated  liquid  viscosity  data  for  the  hydrocarbons  are  primarily 
obtained  from  a  single  source  (69)  and  the  accuracy  of  these  data  are 
not  made  known.   The  high  pressure  viscosity  data  are  less  accurate. 
Thus  for  nitrogen  at  a  pressure  of  250  atmospheres  the  maximum  dis- 
agreement between  workers  is  as  much  as  14%,  and  for  methane  and  ethane 
discrepancies  are  over  10%. 

Saturated  liquid  thermal  conductivities  appear  to  be  of  quite 
good  accuracy  for  most  heavy  fluids,  being  of  the  order  2%.   However, 
for  inert  gases,  for  which  extensive  measurements  have  been  made, 
discrepancies  are  found  to  be  as  much  as  7%.   Of  all  the  three  trans- 
port coefficients,  self-dif fusivity  measurements  are  the  least  accurate. 
Thus  the  methane  data  reported  by  different  workers  show  discrepancies 
as  large  as  18%,  and  for  ethane  differences  are  20%. 
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TABLE  6.2 
SOURCES  OF  TRANSPORT  PROPERTY  DATA 


Molecule 


CH, 


CD, 


Ne 


CO 


^2«6 


Reference 

52,64,68,69, 
115,116,117, 
118,119 

64 

61,62,90, 
107 

64,67,117, 
120 

49,61,64, 
67,68,98, 
108,116,121, 
122 

69,119,123, 
124 

59,64,67,68, 

70,104,108, 

120,121,125 

66 


X 

Reference 

17,27,66,138, 
139,140,141 


66,103,111 
27,66,138 


27,66,100, 
101,102,142, 
143 , 144 


140,145 
66,146 


Reference 


18,58,72 


71 


149,150 


Ar 

15,48,60,61, 

17,27,60,66, 

18 

64,65,67,97, 

100,101,102, 

98,99,104, 

109 

105,108 

co^ 

53,104,126 

138,146,139, 

151,152,153 

147 

CoHc 

52,69,118, 

3  8 

119,124,127, 
128,129 

^-^4^10 

69,130 

148 

108 


TABLE   6.2    (Continued) 


-n 

X 

D 

Molecule 

Reference 

Reference 

Reference 

i-C4%o 

131 

°-S«io 

69 

--^^2 

69,130 

154,155 

C(CH3)^ 

128,130,132, 

156,157 

133 

5   0 

55,69,134 

27,136,147 

156,158,159, 

135 

160,161,162 

Kr 

33,64,65 

17,60,100 

18 

0-CgH^2 

69 

29,147 

156,161,163 

-^6«14 

69,130 

n-C^H^g 

69,130 

^-^8-^18 

69,130 

n-CgH^o 

69 

Xe 

64,65,106 

17,60,100 

18,57 

CCl, 

4 

29,55,128, 

27,136,146 

163,164,165 

134,136,137 
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Figure  6.5  is  a  reduced  plot  of  saturated  liqiid  viscosity  as 
a  function  of  temperature  using  the  improved  correspondence  principle, 
where  the  viscosity  and  temperature  are  now  reduced  according  to 
equations  (o.49)  and  (6.41)  respectively.   For  clarity,  only  representa- 
tive molecules  of  different  molecular  structures  and  weights  are 
included  in  this  figure.   However,  similar  agreement  is  found  for  the 
other  molecules  listed  in  Table  6.1,   The  molecules  used  to  test  the 
modified  correspondence  principle  include  long-chain  molecules  such  as 
nonane,  spherical  molecules  such  as  neopentane,  flat  molecules  such 
as  benzene,  very  light  molecules  such  as  methane  and  neon,  heavy 
molecules  such  as  carbon  tetrachloride  and  isotopic  molecules  such  as 
tetradeuterome thane . 

Figure  6.6  is  a  similar  plot  for  reduced  saturated  liquid 
thermal  conductivity.   The  thermal  conductivity  is  reduced  according 
to  equation  (6.51).   Good  agreement  is  found  except  for  the  data  of 
carbon  tetrachloride,  benzene  and  neon.   The  deviation  for  neon  may  be 
due  to  quantum  effects  as  has  been  pointed  out  in  the  previous  chapter. 
The  poor  correlation  for  carbon  tetrachloride  and  benzene  may  be  due 
to  the  Eucken  effect.   This  effect  appears  to  be  more  important  in  the 
gaseous  state  than  for  liquids,  and  is  discussed  in  more  detail  later. 

Figure  6,7  is  a  plot  of  reduced  self-diffusion  coefficient  for 
saturated  liquids.   The  correlation  is  seen  to  be  poorer  than  for  the 
other  two  transport  coefficients.   This  may  be  due  to  experimental 
inaccuracies. 

Figure  6.8  shows  a  test  of  the  principle  for  viscosity  over  a 
very  wide  range  of  temperatures  and  pressures.   Pressures  up  to  almost 
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1000  atmospheres  are  included.   With  the  modified  principle  good 

~vV      1 
correlation  is  observed  for  the  gas  and  liquid  viscosities  at  P   =0. 

However,  the  correlation  is  poorer  at  higher  reduced  pressures.   This 

may  be  due  to  the  fact  that  c  is  slightly  density  dependent.   Within 

the  normal  liquid  density  range  the  factor  c  is  found  to  be  approximately 

constant. 

The  simple  correspondence  principle  does  not  apply  to  the 

thermal  conductivity  of  polyatomic  fluids  in  the  gas  phase,  because 

of  the  Eucken  effect.   The  Eucken  effect  can  be  largely  corrected  for 

by  plotting  the  gas  thermal  conductivity  as  logW?^"  vs.  T""  ,  where 

o 

X  "  is  arbitrarily  chosen  as  the  value  of  X"   at  T"  =  0.1.   Figure  6.9 
o 

shows  such  a  plot  for  a  nuiriier  of  monatomic  and  polyatomic  molecules. 
Using  this  method  of  presentation,  corresponding  states  correlations 
are  found  to  be  greatly  improved.   It  seems  that  the  Eucken  effect  is 
more  important  in  the  gas  than  in  the  liquid  phase,  for  good  correla- 
tion is  observed  in  the  liquid  phase  (with  the  exception  of  benzene 
and  carbon  tetrachloride)  without  correction  for  the  effect.   A  test 
of  the  modified  principle  for  the  thermal  conductivity  of  polyatomic 
molecules  at  P  =  0.0 (liquid  phase)  and  ?"  =  1.0  is  shown  in  Figure  6.10. 

Tables  6.3,  6.4  and  6.5  give  the  coefficients  of  the  equations 
(5.42),  (5.43)  and  (5.44)  fitted  for  the  curves  shown  in  Figures  6.5, 
6.6  and  6.7  respectively.   Estimated  average  percent  errors  for  these 

equations  are  also  given  in  the  tables.   Table  6.6  gives  the  coeffi- 

1  ~v.- 

The  reduced  pressure  P  used  for  polyatomic  molecules  differs 

from  that  used  in  Chapter  5,  and  is  defined  by  equation  (6.62).   It 

involves  the  term  c.   This  definition  arises  for  the  equation  of 

state,  and  is  discussed  in  the  next  section. 
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cients  of  equations  for  the  curves  in  Figure  6.8  and  Tables  5.7  and 
6.8  give  coefficients  of  the  thermal  conductivity  equations  for  the 
curves  in  Figures  6.9  and  6.10  respectively.   A  list  of  the  values  of 
X"    for  molecules  considered  in  Figure  6.9  is  also  included  in  Table 
6.7. 

6 .4 Comparison  with  Theory  for  Thermodynamic  Properties- 

Recently,  Hermsen  and  Prausnitz  (166)  proposed  a  corresponding 
states  treatment  for  the  thermodynamic  properties  of  liquid  hydrocar- 
bons.  Their  approach  was  based  on  Prigogine's  corresponding  states 
theory  for  r-mers  (50)  .   Th^  defined  the  reduced  temperature  as 


T  =  ,  )\  (6.52) 

and  the  reduced  molal  volume  as 

V  =  ^  (6.53) 

k 

where  V"  is  an  arbitrary  characteristic  volume.   Potential  parameters 

k 

determined  from  second  virial  coefficient  data  as  given  by  Hirschfelder, 

Curtiss  and  Bird  (20)  were  used.   By  fitting  experimental  molal  volume 

data  using  a  high  speed  computer  they  obtained  (s/c)   and  V"  for  a 

k      k 

large  number  of  hydrocarbons.   It  is  of  interest  to  compare  their 

values  of  (s/c)   obtained  from  molal  volume  data  with  the  values  of  c 
k 

obtained  here  from  viscosity  data.   These  quantities,  together  with 
the  potential  parameters  used  in  each  case, are  compared  in  Table  6.9. 
Hermsen  and  Prausnitz  arbitrarily  assumed  the  value  of  (s/c)   for 
methane  co  be  unity.   After  multiplying  their  values  of  (s/c)^  by  a 
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constnnt  factor  chosen  such  that  the  value  of  (s/lO^^  for  methane 
agrees  with  the  value  of  c  for  methane  used  in  this  work,  it  is  seen 
from  Table  6.9  that  the  txvro  parameters  are  approximately  the  same  for 
most  hydrocarbons.   This  suggests  that  the  correspondence  principle 
should  apply  quite  well  to  both  transport  and  thermodynamic  properties 
in  the  liquid  phase  using  a  single  set  of  c  values.   However,  it  may 
be  necessary  to  use  different  sets  of  potential  parameters  in  the  t\Jo 

cases. 

It  is  of  interest  to  examine  the  modified  correspondence 
principle  for  thermodynamic  properties  of  polyatomic  fluids.   From 
statistical  thermodynamics,  the  pressure  is  given  in  terms  of  the 
canonical  partition  function  as 


where  the  partition  function  may  be  separated  into  a  density- independent 

Dart  0  and  a  part  Q   which  depends  on  density, 
^    ^i       ^  CM 

Q=i7QQ  (6.55) 


where 


400   V 
-co   o 

3N   2 

r— ,   CO. 

\  -r^  +   cU^r r.J  (6-57) 


and 
therefore 


CM    Z    ^™       -i    1^ 

i=l 

U(r,...r^.)-.  V  *(^)  (6.58) 
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P  =  kT 


V    V 
^(In  /. . .  /  e  ^     dr^. . .dv^) 


bv 


(6.59) 


N  cpf 

-K»      -foo        y   ik 


Since 


1=1  2mkT 
e  cl£i''''^£»T   is  independent  of  density. 


-00    -00 


Equation  (6.59)  may  be  written  in  terms  of  reduced  variables  to  give 


ce  ~* 


V    V 


S(ln 


exp 


o   o 


L   T 


0(rt.) 


i<J 


hv* 


or 


where 


P  =  P  (T  ,V  ) 
ce 

CT 


T 

(6.60) 
(6.61) 


(6.62) 
(6.63) 


Pitzer  (94)  has  pointed  out  that  the  simple  2-parameter 
corresponding  states  law  fails  most  noticeably  for  the  liquid-vapor 
region,  and  suggested  using  experimental  vapor  pressure  data  to  obtain 
a  third  parameter.   Figure  6.11  shows  a  plot  of  reduced  vapor  pressure 
P  (T  )  using  the  simple  2-parameter  corresponding  states  principle; 
substantial  deviations  are  seen  to  occur,  and  these  are  ascribed  by 
Pitzer  to  acentric  interactions.   Figure  6.12  shows  a  similar  plot  of 
P  (T  )  using  the  new  correspondence  principle,  with  the  values  of  c 
obtained  from  liquid  viscosity  data.   Potential  parameters  used  were 
smoothed  values  from  gas  viscosity  data.   The  correlation  is  considerably 
improved.   Further  improvement  can,  of  course,  be  obtained  if  c  and 
potential  parameters  are  derived  from  thermodynamic  data. 
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CHAPTER  7 
CONCLUSIONS 

The  free  volxur.e  equation  provides  a  useful  basis  for  pre- 
dicting viscosities  and  self-diffusion  coefficients  of  liquids  at 
high  densities,  where  most  other  correlations  are  unsuccessful.   The 
fact  that  both  V  and  E  values  vary  in  a  corresponding  states  way 
indicates  that  the  equation's  success  does  not  arise  only  from  the 
availability  of  three  adjustable  parameters.   The  simple  expressions 
given  for  A  by  Eyring  and  others,  however,  do  not  seem  to  be  general- 
ly applicable.   The  free  volume  theory  seems  to  provide  somewhat  more 
accurate  values  of  transport  properties  than  does  the  direct  corres- 
ponding states  approach,  but  has  several  limitations.   In  particular 
it  applies  only  ac  high  density  (p>  2p  ),  and  cannot  be  used  to  obtain 
thermal  conductivities. 

The  direct  corresponding  states  approach  is  of  more  general 
application,  and  in  the  form  presented  in  Chapter  6  it  may  be  used  for 
all  transport  properties  and  over  a  wide  range  of  densities.   An 
attractive  feature  of  the  principle  presented  for  polyatomic  fluids 
is  that,  although  a  third  parameter  is  introduced,  the  transport 
coefficients  are  expressed  as  functions  of  only  two  independent 
variables  (T  and  P  ) .   In  view  of  the  simplicity  of  the  proposed 
model,  it  may  at  first  sight  seem  somewhat  surprising  that  the 
correlation  is  so  successful.   It  should  be  remembered,  however,  that 
the  model  is  only  used  to  obtain  the  form  of  the  reduced  relation- 
ships, and  not  to  directly  calculate  transport  coefficients.   In  this 
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forra  the  correspondence  principle  applies  to  a  wide  variety  of 
substances,  including  fluids  composed  of  relatively  complex  molecules 
X'/hich  are  of  interest  to  engineers.   It  is  easy  to  apply,  the  only 
parameters  required  for  prediction  being  T  ,  P  ,  M,  CD  and  c.   The  first 
four  parameters  are  available  in  the  literature,  and  the  latter  may 
be  obtained  from  one  experimental  measurement  of  any  one  of  the 
transport  properties.   I^Jhen  no  such  measurements  are  available  it 
may  be  evaluated  using  thermodynamic  data  (e.g.,  vapor  pressure). 

In  the  light  of  the  promising  results  of  this  investigation, 
it  is  felt  that  the  free  volume  theory  and  the  corresponding 
states  principle  are  very  useful  methods  for  the  prediction  of  phy- 
sical properties  of  dense  and  condensed  fluids.   For  future  work 
the  following  studies  are  recommended. 
7  .1 Free  Volume  Theory 

7.1.1  Extension  to  Mixtures 

In  the  course  of  this  work  the  free  volume  theory  has  been 
tested  for  the  viscosity  of  binary  mixtures  of  simple  molecules 
(167) .   For  the  mixtures  tested  it  is  found  that  the  predicted 
viscosity  is  within  the  experimental  accuracy  over  the  entire  range 
of  compositions.   In  the  prediction  of  the  mixture  viscosity  no 
mixture  data  are  required.   For  future  studies,  the  free  volume 
theory  may  be  extended  to  binary  diffusion  coefficients  and  viscosity 
of  multicomponent  systems. 

7.1.2  Electrolyte  Solutions 

Podolsky  (168)  has  utilized  the  Eyring  cell  theory  to  pre- 
dict the  self-diffusion  coefficient  of  electrolyte  solutions. 
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However,  no  attempt  has  yet  been  made  to  apply  the  free  volume 
theory  to  the  transport  properties  of  electrolytes.   Based  on  the 
promising  results  demonstrated  by  the  free  volume  theory  for  the 
viscosity  and  self-dif fusivity  of  pure  fluids,  it  is  believed  that 
this  theory  may  serve  as  an  improved  basis  for  predicting  transport 
properties  of  electrolytes. 
7.2 Corresponding  States  Principle 

7.2.1  Mixtures 

The  corresponding  states  principle  may  be  derived  in  the  same 
manner  as  for  pure  substances,  by  starting  with  the  time  correlation 
function  equations  for  the  transport  properties  of  mixtures  (169-171). 

7.2.2  Polar  Substances 

In  this  work  only  nor.polar  substances  are  considered.   It  will 
be  useful  to  extend  the  correspondence  principle  to  polar  substances. 
This  can  be  done  by  introducing  suitable  reduced  dipole  or  quadrupole 
moments , 

7.2.3  Fused  Salts 

Transport  properties  of  fused  salts  are  very  important  for  a 
number  of  practical  applications.   As  measurements  of  transport  pro- 
perties of  fused  salts  are  difficult,  it  will  be  of  interest  to  extend 
the  corresponding  states  principle  to  these  substances.   Such  liquids 
are  composed  of  ions  which  may  often  be  regarded  as  roughly  spherical 
charged  particles. 

7.2.4  Thermodynamic  Properties 

A  further  extension  of  the  proposed  principle  for  polyatomic 
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molecules  will  be  a  detailed  study  of  its  application  to  thermodynamic 
properties  in  the  der.Sj  i/iu  condensed  phase.   A  preliminary  study  shows 
that  the  proposed  correspondence  principle  greatly  iu.proves  the  correla- 
tion for  vapor  pressure.   With  suitable  reduced  expressions,  it 
seems  likely  that  the  correspondence  principle  may  apply  to  all  other 
thermodynamic  properties. 
7.2.5   Prediction  of  c 

Although  an  approximate  expression  has  been  obtained  for  the 
factor  c  from  the  derivation  of  the  correspondence  principle  for 
polyatomic  molecules,  no  extensive  study  of  the  behavior  of  this 
parameter  has  been  made.   The  parameter  c  is  a  function  of  the  moment 
of  inertia,  mass  and  potential  parameters  of  the  molecule,  as  well 
as  the  density  and  temperature.   A  useful  extension  of  this  work  will 
be  the  study  of  the  properties  of  the  parameter  c  and  methods  of 
predicting  it.   Any  more  detailed  examination  should  include  a  study 
of  the  density  dependence  of  this  parameter. 
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APPENDIX  1 

CHUNG'S  DERIVATION  OF  THE 
FREE  VOLUXE  EQUATIONS 

Chung  (46)  considered  a  aiolecule  to  spend  a  large  fraction 
of  the  time  oscillating  about  an  equilibrium  position  in  a  cell. 
The  m.olecule  can  leave  the  cell  to  jump  into  a  neighboring  cell 
only  when  the  following  b<!0   events  occur  simultaneously:   (1)  a 
free  volume  greater  than  the  minimiom  value  v  ,  which  can  accomodate 
the  molecule  is  available,  and  (2)  the  molecule  has  an  activation 
energy  greater  than  the  minimimi  energy  barrier  e   separating  two 
equilibrium  positions. 

Consider  a  system  of  N  molecules.   The  total  energy  is 

considered  to  be  made  up  of  energy  levels  e-,,  e  .  ...e.,  such  that 

N  ■'  molecules  are  in  region  e  ,  N  ^    in  e  ....  and  so  on.   In 

addition  it  is  assumed  that  a  range  of  values  of  free  voliome  v  , 

V  . . .V   is  available  to  each  molecule.   Thus  the  distribution  of 
2    k 

molecules  among  the  possible  values  of  free  volume  and  energy  may 

be  represented  schematically  by 

Energy  level       e^  e„  e. 

12  J 

Free  volume  v,v...v  vv...v  vv...v 

12         k  12         k  12k 

Number  of       „(1)„(1)    .,(1)   .,(2)„(2)     (2)    ^0)^0)   ^,0) 
molecules       N   N,   .  .  .N      N^   N   .  .  .N       N   N   ...K 
1^      k      12      k       12      k 

Thus  if  W  represents  the  number  of  distinguishable  ways  of  assigning 

N  molecules  to  the  different  energy  levels  such  that  N^  '^  molecules 

have  energy  e  ,  N^  "^  have  energy  e^  ...  and  so  on,  W  may  be  written 
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as 


\-   1^     ^(i),  (Al.l) 

i=l 
Similarly,  if  W  is  the  number  of  ways  of  assigning  N    molecules 

in  V  ,  N  "^  in  V   ...  etc.,  where  i  =  1  to  j ,  it  follows  that 
12       2 


^    ^    ,       N^^^:     1  r      N^J^' 


p=l  p  p=l  p 

1=1 


j    k  ^ri),  (A1.2) 

i=l  p=l  P 
The  total  number  of  configuration  becomes 


W  =  WW 


iN(^>: 


n:         i=i 


TT  N^  '  !    TT    7T  N^  '  I 

i=l       i=l  p=l  P 


n: 


"  j    k^,(i),  (A1.3) 

i=l  p=l  P 
Making  use  of  Stirling's  approximation  we  may  write 

j   k 


InW  =  NlnN  -   )   V   N^^^^lnN^^^ 


(A1.4) 


i=l  p=l 
The  statistical  mechanical  formulation  of  entropy  S  is  given  by 


S  =  klnW      .   , 

=  k[NlnN  -  V  Vw^^-'lnN^^'*]  (A1.5) 


_  -^  P     P 

i=l  p=l 
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The  system  is  isolated  and  is  considered  to  be  at  thermodynamic 
equilibriiim.   Thus  the  entropy  is  a  maximum,  and 


dS  =  0 


j   k     ^  j   k 

=  k)  V   — 7TT  I  NlnN  - 

1=1  p=l   p  m=l  n=l 


N^^^nN^^^ 
n     n 


dN 


(i) 


J   t^ 


=  -  k 


?^ 


i=l  p=l  P 


P     P 


dN 


.(i) 


J   k 


-k  yy  [InN^^)  -M]d4^^ 


£_  ^        p 

i=l  p=l 


(A1.6) 


or 


j   k 
y  y   [inN^^^  +  l]dN<^^^  =  0 

/    /  D  D 


z_ 

i=l  p=l 


(A1.7) 


The  system  is  subject  to  the  restraints  of  constant  number  of 
molecules,  constant  total  free  volume  and  constant  energy,  so  that 


j   k 


or 


P 
i=l  p=l 

1 1  ^A''  =  ^ 


i=l  p=l 

j     k 

i=l  p=l 

^f 

j      k 

Z.  L     p 

i=l  p=l 

0 

j      k 

L.  L     ^   p 

i=l  p=l 

=  0 

j      k 

y  y  V  dN(^>  =  0 

L  L     p   p 

i=l  p=l 


(Al.S) 


(A1.9) 


(ALIO) 


(Al.ll) 


(A1.12) 


(A1.13) 
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Using  the  Lagrange  undetermined  multiplier  procedure  we  may  multiply 
equations (Al. 11),  (A1.12)  and  (A1.13)  by  (a-  1),  /3  and  A  respectively, 
and  add  the  resulting  equations  to  equation  (A1.7),  to  give 
J    k 

^    >     [a  +  InN^^-*  +  /Se.  +  Av  ]  dN  ^^^    =   0 
L.       i.  p     ^  1     p'    p 

i=l   p=l 

Choosing  suitable  values  for  a,  jS  and  A,  one  may  write 


(A1.14) 


or 


a  +  InN^^-*  +  fie.  -rAv   =0  (A1.15) 

pip. 

-.^x    -(a  +  pe.  +  Av  ) 
Np   =  e        ^     P  for  all  p,i     (A1.1&) 

The  probability  distribution  of  free  volume  and  energy  may  be 
written  as 

(l)       D 

Pp     =  -— 

-(a  +  fie .  +  Av  ) 
1     p' 
e  ^ 


Y"!  r-i   -  (a  +  /3e ,  +  Av  ) 

^   )   e        1     P  (A1.17) 

i=l  p"^l 

From  equation  (A1.16),  on  summing  over  p  and  i,  we  have 


X     X       -  (a  +  /5e .  +  Av  ) 
N  =  V  V  e      ^     ? 


i=l  p=l 
Solving   for  a  gives 


(Al.lS) 


-Q.  N 

e 


^  V    -<^^i-^%) 


e 


(A1.19) 


1=1  p=i 
Making  use  of  equations  (ALIO),  (A1.16)  and  (A1.19),  and  allowing 

the  values  of  free  volume  v^^  to  tend  to  the  continuum  limit,  we  have 
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X  ^     -(«e.  +Av^) 

1       D 

P 

N  )   e      )   V  e   '^ 

_   i=l D=2 

=  j    ^    k   ,  (A1.20) 


I  "^  I 


e   P 


thus  i='l      i=l 

V     J  ve   dv 


f   N       r=^  (A1.21) 


a'^ 


-Av  ~ 
e   dv 

where   v  =  fluctuating  free  voliome 

or  A=—  (A1.22) 

^f 

Similarly  combining  equations  (A1.9),  (A1.15)  and  (A1.19)  gives 

^  -Be. 

^   1 


V 


e.e 

1 


^    X  (A1.23) 


-Be. 
e        "- 


i=l 
:t  is  obvious  from  equation  (A1.23)  that  /3  has  its  usual  significance. 


that  is 

'^  "  kT 


^=h  <A1.24) 


Equ:  tion  (A1.17)  can  nox^  be  written  as 
P^'^  =  P(vp)p(e.) 

e:-:p(-v  /v^)ex?(-e, /kT) 


>  V  exp(-v  /v  )exp(-e./kT) 

Z_;  4_        P   ^       ^  (A1.25) 

1=1  p=l 


142 


Since    the  values   of  v      are   assumed   to  be   continuous,    the   sicnmation 

P 

can  be  replaced  by  an  integral  and  the  subscript  p  dropped,  so  that 


p(v)  =-expf-  -J 


(A1.26) 


To  get  an  explicit  expression  for  the  energy  distribution,  Chung 
assumed  the  molecules  to  be  harmonic  oscillators.   In  the  classical 
limit  (172) 


?(e)  = 


1  \   /'e_ 

\kT 


2kT 


exp 


e 
kT 


(A1.27) 


The  probability  that  a  molecule  simultaneously  has  an  energy  greater 
than  some  value  e  and  a  free  volume  greater  than  some  value  v   is: 


P(v  >  v  ,e  >  e  )  = 

O  V 


p(v)p(e)dedv 


v   e 
o   V 


,  ,  e  \  2 

1  /   V 


2  \.  kT 


+ 


V 

kT 


+  1 


exp 


r  /e    V 
\   kT   v^ 


(A1.28) 


Based  on  the  assijmption  that  self-dif fusivity  and  fluidity  are 
directly  proportional  to  this  probability,  one  may  write  for  the 
expressions  of  viscosity  and  self-dif fusivity 

E     V 


r  1 


(A1.29) 


and 


D  =  Bexp 


/  E      V 

f   _Z  ^ 2_ 

V  RT   V  -  V 


(A1.30) 


where 
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V  .  =  V  -  V 

t        o 
V  =  Nv 


V  =  Nv 
o     o 

E  =  Ne 

V       V 


and  A  and  B  are  assumed  roughly  constant.   The  form  of  the  pre- 
exponential  factors  is  not  made  clear  in  the  derivation.   The  pre- 
exponential  term  appearing  in  equation  (A1.2S)  depends  on  the 
assumed  form  of  the  molecular  energy.   Thus  if  the  molecules  are 
assumed  to  be  in  a  potential  well,  rather  than  behaving  as  oscil- 
lators, a  different  pre-exponential  factor  is  obtained. 


APPENDIX   2 
FURTHER   TEST  OF   PROPOSED  FREE  VOLUME   THEORY 
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o Experiment    (49) 
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Figure  A2.1.   Test  of  Free  Volume  Theory  for  Liquid  Nitrogen 
Viscosity. 
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APPENDIX  3 

SOLUTION  OF  EQUATIONS  OF  MOTION 
FOR  OSCILLATOR 


From  equations  (6.15)  and  (6.17)  we  have 


m 


d^x 


dt 


2-^^l^  +  ^2^=° 


(A3.1) 


and 


2 

d  -^ 


d-c 


z     ZZ  X         IZ 


(A3. 2) 


In  terms  of  differential  operators  one  may  write 


and 


(mD  +  k   )x  +  k.„^  =  0 
11      12  X 


(ID  +  k.,)?  +  k.^x  =  0 
ZZ     X    12 


(A3. 3) 
(A3. 4) 


Eliminating  Y.from  equations  (A3. 3)  and  (A3. 4)  gives 


(ID^  +  k22)(mD^  +  k^^x  -  k^^^  =  ^ 


(A3. 5) 


Expanding  equation  (A3. 5)  we  have 


[imD^  +  (k^^m  +  k^^I)D^  +  k^^k^,  -  k^^^^  =  ° 


(A3. 6) 


Equation  (A3. 6)  may  be  solved  to  give  roots  of  D  as 


n 


(k22in+k^,^I)+y(k22m+k^^I)^  -  4lm(k^^k22-k^,) 


1/2 


(A3. 7) 


Let 
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^1  = 


(k^.,ni+k^^T')    -\/ (k.^./>ilk^,  l)"^    -    4Im(k^^k^2"'^i9)      " 


1/2 


2Im 


(A3. 8) 


0),    = 


(k2om+k^^I)   +/(k22m+k^^I)^   -   4lm(k^^k22-k^2) 


1/2 


2Im 


(A3. 9) 


Therefore 


-ico  t  +ico  t  -iCJ„t  ico^t 

X  =  A.e  +  A„e  +  A„e  '  +  A,  e 

12  3  4 


(A3. 10) 


Similarly,    eliminating  x   from  equations    (A3. 3)    and    (A3. 4)    gives 


(mD^   +  k-^)(ID^   +  k„)^  -   k^.Y  =   0 
li  2.1      yi         1/    X 


(A3. 11) 


Equation  (A3. 11)  is  expanded  to  give 


[iraD^  -r  (k^2™  +  k^^I)D^  +  k^,^k22  -  k^2^^  "  ° 


(A3. 12) 


Solving   for   the   roots   of   D  we  have 


D  =    '    i 


-    /  ~  2/        1/ 

j-    (k22m+k^I)    +  ^/(k2^m+k^I)      -    4Im(k     k^^-k^)      -, 


2Im 


(A3. 13) 


Making  use   of  equations    (A3. 13),    (A3. 8)    and    (A3. 9), 


-iaj  t  ico,  t  -iCD^t  i6D2t 

H'     =  E,e  +  B.e  +  B„e  +  B,  e 

X  1  2  3  4 


(A3. 14) 


Since  v/e  have  assumed  that  the  translational  and  angular  oscillations 
represented  by  equations  (A3.1)  and  (A3. 2)  are  harmonic  and  that  the 
oscillations  are  in  phase,  then  60  ,  which  represents  the  angular 
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velocity  of  a  harmonic  motion,  must  be  equal  to  o)^  •   This  is  true 

if 


(k22m+k^^I)^  -  4lm(k^^k22-k^2^  =  °  (A3. 15) 

Equation  (A3. 10)  now  reduces  to 

X  =  A'e   ^  +  A^e""'^!*'  (A3. 16) 

-icd  t 
Since  e  =  cosOO  t   -    isimo  t  (A3. 17) 

i(X>- 1 
and  e  =  cosCD  c   +  isinO).  t  (A3. 18) 

1  ->• 

Equation  (A3. 16)  may  be  written  as 


X  =  A'cosO),  t  +  A"sinJO,  t  (A3. 19) 

1  i. 


where 


A'  =  A|  +  A' 
A"  =  i(A'  -  Ap 
Similarly  equation  (A3. 14)  may  be  written  as 


^  =  B'cosco.t  +  B"sinw,t  (A3. 20) 

X         1  i 


APPENDIX  4 
HAMILTONIAN  FOR  THE    TIIREE-DIMENSIONAL  OSCILLATOR 


X 


Figure  A4.1.   A  Three-Dimensional  Oscillator  Model 

Consider  a  molecule  oscillating  in  an  isotropic  harmonic 
manner  in  a  3-dimensicnal  cell.   At  any  particular  position  _r  of  the 
molecule,  it  has  a  linear  velocity  V(_r)  and  an  angular  velocity 
co(_r,_0) ,  where  0  represents  the  Eulerian  angles. 

One  may  resolve  _r  and  V(_r)  along  the  x,  y  and  z  axes  as 


r  =  xe  +  ye  +  ze 


(A4.1) 


V(r)  =  V^(x)e^  +  V  (y)e2  +  V^(z)e3        (A4.2) 


Let  to(_r,_9)  be  made  up  of  two  parts;  a  constant  angular  velocity 
CJ°  (_r  =  O,_0)  together  with  a  harmonic  angular  oscillator  which  is 
center  of  mass  coordinates  dependent,  a)(_r)  ,  superimposed  on  it.   In 

other  v7ords  we  assume 


(^(r,e)   =  a.(r)  +  (a   ( r  =  0,0) 


(A4.3) 
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Resolving  this  along  the  3  Cartesian  coordinates  we  have 

a^(r,_§)  =  a).(x)e  +  CD  (y)e  +  O^Jz)e     +   a)°(^)e 
X    X    y    z    ^    J        J. 

+  a3°(0)e  +6D°(9)e  (A4.4) 

y  — 2   ^  ~3 

The  angular  velocities  at  the  cell  center  can  be  expressed  in  terms 
of  Eulerian  angles  as  (173) 

CD°(6;)  =  0COS0  +ijj   sin0sin0  (A4.5) 

X 

CO°(e)   =    0sin0  -  •^sin9sin0  (A4.&) 

GO°(0)   =  0  +  7//cose  (A4.7) 

z  ~" 

Suppose  that  the  niolecule  is  oscillating  in  a  particular 
direction  with  instantaneous  velocity  V  at  any  position  r,  and  at 
the  same  time  it  has  a  harmonic  angular  oscillation  which  is  in 
phase  with  the  linear  one. 

We  now  consider  the  components  of  V  and  Cd(r,6)  .     We  pair 

CV  CO  ) ,  (V  ,03  )  and  (V  ,co  ) ,  and  consider  them  as  the  linear  and 
X  z     y  X        z  y 

angular  velocities  of  3  oscillators,  one  along  the  x-axis,  one  along 

the  y-axis  and  one  along  the  z-axis. 

These  oscillators  may  be  considered  separately.   The  oscillation 

associated  with  the  x-axis  will  be  as  shown  in  Figure  6.3.  We  define 

the  angles  associated  with  co  as 
°  X 

t 

e  (t)  =    r  CO  (x,0)dt 

■^-      J        ^  (A4.8) 

t 
o 
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^   (t)  =  /  co,(e)dt  (A4.9) 


t 
o 

t 


=  r  o)  (x^ 


^  (x)  =   /  O)  (x)dt  (A4.10) 

X  ,'     X 


t 
o 


The  potential  and  kinetic  energies  are  expressed  as  before  as 


1   .2    1   '2 
K  =  ^  mx  +  ^  I^     •  (A4.12) 

X   2       2   X 

The  equations  of  motion  and  their  solution  are  nov/  identical  to  those 
given  in  section  6.2.   The  final  equations  for  the  potential  and 
kinetic  energies  associated  with  the  x-direction  are 

=  I  c[k^^  +  Hk^^'^^  (A4.13) 


K  =  ^  cxax-  (A4.14) 

where  H  and  c  are  defined  by  equations  (6.32)  and  (6.37).   Similarly 
for  the  y  and  z   axes  we  have 

^/y.\)    =  ^   c[k^^  +  Hk^2^y2  (A4.15) 

U^(.,Y^)  4  c[k^^  +  Hk^2^z2  (A4.16) 

K  =  -J  cmy^  (A4.17) 
y   2   ^ 
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K  =  —  cmz 
z   2 


(A4.18) 


The  Hamiltonian  for  the  center  of  mass  coordinates  becomes 


H 


CM 


N 

)     K.   +  K.   +  K.   +  U.   +  U.   +  U.      +  U(0,0) 
^_,     IX    ly     12     IX    ly     iz       v  j^-' 

i=l 
N 

i=l  ~  -^ 


=  Vfc 


2,2,2 


i=l 


zm 


(p.   +  p.   +  p.  )  +  cU(r r  ,) 


(A4.19) 
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